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ABSTRACT 



A new consideration in the growth mechanism of whiskers is exam- 
ined. The heat flow due to the shape of the crystal and the thermo- 
dynamics of the growth reaction is believed to be an important factor 
in the unidirectional growth observed in whiskers. A computer pro- 
gram was formulated to model the heat flow as well as a method to 
calculate the time necessary for the temperature changes to occur. 

It was found that there is substantial reason to believe the whisker 
tip is significantly cooler than the sides. This phenomenon is referred 
to as the cold tip mechanism which provides a preferential growth site 
for reaction at the tip. 
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I. INTRODUCTION 



The requirements of today's technology have brought about the de- 
mand for more high strength, light-weight materials. In the past few 
years, fiber reinforced composite materials have shown promise in 
fulfilling this demand. Of concern in this paper is the growth of a class 
of fibers called whiskers. The terms "fiber" and "whisker" have been 
defined by Rauch, Sutton, and McCreight in Ref. 1 as: 

Fiber - Any material in an elongated form such that it 
has a minimum length to a maximum average transverse 
dimension of 10:1, a maximum cross-sectional area of 
7. 9x10" 5 in .2 and a transverse dimension of 0.010 in. 

Whisker - Any material that fits the definition of a fiber 
and is a single crystal. 

Whiskers are exceptionally suited to this application due to their high 
strength. This strength has been shown by Brenner (Ref. 2) to app- 
roach the theoretical cohesive strength of solid matter, which lies be- 
tween three and seventeen per cent of the modulus of elasticity. This 
remarkable property has been attributed to the crystalline perfection 
and small dimensions that minimize defects or dislocations which 
weaken bulk materials. Frank has estimated that a bulk material of 
comparable size to a whisker would contain 10® dislocations (Ref. 3). 

Currently, the major reinforcing whisker in use is of die ceramic 
type. Metallic whiskers have not been used for engineering purposes 
to as great an extent due to problems of fabricating them into a matrix. 
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In general, a metallic whisker undergoes metallurgical reaction with 
the matrix at high temperature and during fabrication. Dislocations 
are easily introduced into most metallic whiskers, and, in comparison 
to ceramic whiskers, possess a low modulus and high density. 

To solve the problem of fabrication and production, an understanding 
of the basic mechanism for metallic whisker growth is necessary. 

There are two generally accepted mechanisms for whisker growth. 

The first mechanism proposed is based upon the screw dislocation 
theory (Ref. 3). Basically, this mechanism postulates a permanent 
crystal imperfection at the tip of the whisker. The dislocation be- 
comes a preferred growth site for the material which diffuses up the 
sides of the whisker. The second common theory is that of the vapor- 
liquid- solid (VLS) mechanism (Ref. 1). Here an impurity becomes a 
site for condensation of the material from the atmosphere. A liquid 
droplet forms, becomes supersaturated, and the growth step occurs 
when the solid precipitates from the liquid. 

Several laboratory experiments were also carried out in the course 
of this investigation concerning growth of whiskers. Growth of iron and 
copper whiskers was examined since these two metals are used widely 
for engineering purposes, and the problem under consideration was the 
separation of defect free whiskers from the substrate. This investiga- 
tion involved growth of whiskers on granular substrates such as char- 
coal, while other experiments delt with transport of the reaction vapor 
from its source to nucleating sites adjacent to the starting materials. 
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In the course of experimentation it was found that the mechanisms pro- 
posed earlier had some shortcomings when used to explain the results. 
Another process was proposed to explain the observed growth in terms 
of heat transfer. 

The new theory was called the cold tip mechanism. This scheme 
proposed that the heat flow in the rod shaped whisker made the tip 
cooler than the sides. The cooler tip provided a position for concen- 
tration of the reactant molecules which had migrated to the whisker 
surface. Consequently, the tip became a preferred growth site. An 
accompanying phenomenon called thermal runaway was predicted. This 
term described the effect of an endothermic reaction withdrawing heat 
from the metal thqs cooling the surface at the tip. Consequently, the 
tip should become an even more preferable growth site. 

To examine this hypothesis, a mathematical model of a cold tip 
whisker was formulated. The finite difference technique was used to 
solve the differential equations describing conduction and convection. 
Results of this treatment were used to show the necessary conditions 
for the cold tip phenomenon to occur. 
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II. TRADITIONAL POSTULATED GROWTH MECHANISMS 



A. SCREW DISLOCATION THEORY 
1. Steps Leading to its Creation 

A theory of crystal growth from the vapor phase was 
initially developed by Volmer in 1922. His early work gave him the 
impression that the mechanism involved an adsorbed layer, and from 
this, the two-dimensional nucleation mechanism was formulated. 
Essentially this theory postulated growth occuring as follows: first, 
consider a flat crystal surface in contact with a vapor. Assume that 
this surface is partially covered by another layer. If the pressure of 
the vapor is raised by an amount called Ap above the equilibrium 
pressure the layer will grow at a speed proportional to Ap until the 
surface is covered. To start a new layer, a two-dimensional nucleus 
(critical nucleus) must be formed. The number of critical nuclei 
created per second is proportional to exp (-Ag/kT), where Aq is the 
activation energy for nucleation. Burton and Cabrera in Ref. 4 cal- 
culated a corrected Aq by: 

0 = LkT Inc* A 0 = 0 2 /(kT lnc<) 

where c< is the saturation ratio, defined as the ratio between the 
actual concentration in the vapor to the equilibrium value, 0-the energy 
of interaction between nearest neighbors, and L the number of 
molecules in the critical nucleus. Using the above equations, 
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reasonable values of Aq were calculated to account for self-nucleation 
of new crystals from the vapor as well as from new liquid drops at 
saturation ratios of 10 and 5 respectively (Ref. 3). However, 1. 5 was 
calculated to be the lowest c< possible to obtain reasonable rates of 
crystal growth by means of two-dimensional nucleation. This result 
deviated markedly from observed crystal growth which occurred in 
environments where super saturation ratios were of the order of 1.01. 

2. Screw Dislocation Mechanism 

To account for this discrepency, F. G. Frank (Ref. 3) 
introduced the screw dislocation mechanism. The screw dislocation 
was pictured as a spiral staircase with the steps made up of crystals 




Fig. 1 Screw dislocation surrounded by growth cone. 



within the lattice at the tip of the whisker. These indentations formed 
preferential growth sites. The dislocation was self-perpetuating in 
that as the steps were filled in by new crystals, and the surface of the 
new growth was of the same orientation as the original growth site. At 
the same time, the screw dislocation would justify the one- dimensional 
growth observed in whiskers. This mechanism has been shown to 
occur in liquid phase crystallization by the observance of growth 
spirals of the above configuration on paraffin. However, conclusive 
experimental evidence that screw dislocations exist in metal whiskers 
grown from the vapor phase has yet to be found. 



B. VAPOR- LIQUID- SOLID (VLS) MECHANISM 

Of the two major whisker mechanisms, VLS has been more recent- 
ly introduced. First, some of the observations which led to the pro- 
posed mechanism for silicon were (Ref. 1): 

1. Silicon whiskers were dislocation free. 

2. Certain impurities were essential for growth. 

3. Whisker growth was a two- stage process 

consisting of : (a) A fast initial extension in 
length (leader growth) followed by (b) a slow 
increase in thickness (layer growth). It was 
also found that impurity concentration in 
leader growth was greater than in the layer growth. 

4. Rapid initial extension in length occurred by 

addition of material at the tip. 

5. A "liquidlike" globule was often observed at the 

tip of the whisker. 

The initiating step in whisker growth by this mechanism is a "seed" 
made up of an impurity which has melted. The liquid surface has a 
higher accomodation coefficient than the side of the whisker and be- 
comes a preferred growth site. The liquid becomes supersaturated 
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with the deposited material from the vapor, and the material precipi- 
tates f or the growth step. 

Since the process is essentially a condensation reaction, the 
growth rate J can be estimated from gas kinetics as: 

J = a <r p 0 (2Trm kT)"2 

Where: 

a = accomodation coefficient - the fraction of impinging 
atoms which deposit on the surface. 

(y = (P^Po) - defined as supersaturation where p is the 
vapor pressure and p 0 is the equilibrium vapor 
pressure of the solid at the temperature T. 
m = mass of molecule, 
k = Boltzmann's constant. 

The accomodation coefficient depends on the perfection of the sur- 
face as well as on 0“ . Thus for a "perfect" crystal no spaces for 
molecules to attach themselves to the whisker exist, and the accomo- 
dation coefficient is small. The rougher the surface the higher the 
coefficient therefore a liquid can be considered the ideally "rough" 
surface since it is composed of accomodation sites only interatomic 
distances apart. Consequently, the VLS mechanism should approach 
the ideal growth rate even at lower supersaturations. 

For VLS to operate, the liquid must be stable, imposing a condition 
upon the diameter of a VLS whisker. The stability of a droplet of 
curvature k in its own vapor depends on c by the expression for the 
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minimum curvature: 



r min = 2 (<rLV)v L 
RT lnCT 



Where: 

( T LV is the liquid - vapor interfacial energy. 

Vl is the liquid molar volume. 

R is the gas constant. 

Also, the VLS mechanism necessitates having an impurity meeting the 
following specifications: (1) It must from a liquid solution with the 
material to be crystallized at the growth temperature. (2) The 
distribution coefficient defined as k = C g / Cp where C g and Cj are the 
solubilities of the solid and liquid respectively, must be less than one. 
(3) The equilibrium vapor pressure of the impurity should be small. 
Evaporation of the impurity lessens the size of the liquid globe, thus 
decreasing the whisker diameter. (4) The seed must be inert to 
chemical reaction with the environment. (5) The agent should have 
the proper wetting characteristics as determined by the interphase 
energies. (6) For controlled unidirectional growth, the solid- liquid 
interface must have a well defined crystal structure. 

The VLS mechanism has much more substantial evidence than the 
screw dislocation mechanism. However, there is no proof that this 
applies to the growth of the pure metallic whiskers. 
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III. EXPERIMENTAL 



A. GROWTH OF IRON WHISKERS 

The method used for growing iron whiskers was introduced by 
Brenner' (Ref. 5) in 1956. The iron was produced in whisker form from 
the reduction of its halide by hydrogen at a temperature above the 
melting point of the salt. 

FeX 2 + H 2 -> Fe + 2 HX 
Where X = a halide. 

The general procedure was as follows: the starting materials were 
prepared and placed upon a steel liner which was used so that the re- 
sultant growth wouldn’t be damaged when removed from the boat as 
well as to keep the boat clean. Liner and boat were placed in the 
cooling section and the rubber stoppers positioned as shown in Fig. 2 
so that the tube was air tight. An inert atmosphere was created by 
flushing the system with helium at a rate of 1. 0 cubic foot per hour 
(CFH) for ten minutes on the first run. For subsequent runs flushing 
time was decreased to approximately five minutes due to the smaller 
amount of oxygen present within the quartz tube. On the initial run of 
a series care was also taken to flush the tubing leading to the hydrogen 
tank to remove any air present. After the chamber was flushed, the 
hydrogen flow was introduced. A short interval of time was allowed to 
pass before the boat was moved into the hot region of the chamber so as 
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Apparatus for growth of whiskers. 



to have a hydrogen atmosphere present for the initial reaction. When 
the desired time for the growth had passed the boat was moved into the 
cooling section, and the system flushed again with helium. The boat 
was removed, and stopper replaced so as to prevent an excessive 
amount of air from entering the chamber. 

This was the general procedure used for growing whiskers. How- 
ever, several modifications were made and will be noted in the following 
experiments. A word about nomenclature in the labeling of experiments 
is appropriate at this time. The convention used in the following pages 
is shown by the example, 2(b)4. 2 stands for the main heading of the 
series; here it stands for the series dealing with "Experiments with 
Ferric Oxide and Ammonium Bromide." (b) serves to differentiate the 
experiments done under this heading. 4 stands for a slight variation of 
the experiment. Usually this entails changing a small amount of 
material or flow rate of gas. Upon seeing the third number, one should 
try to see a trend being formulated within the series. The series often 
involves a trial and error process in an attempt to find an optimun 
effect in accordance with the goal as stated in the purpose under the 
main heading. If no number is present after the letter the experiment 
was unique. 

1. Experiments in a Helium Atmosphere 
Purpose: 

According to Kittaka and Kaneko in reference 6, iron 
whiskers were grown in great quantities by mixing FeCl 2 • 6 I^O or 
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FeBr2 • 2 H 2 Q with carbon black in ratios ranging from 1 to 100 
weight per cent. However, the specific type of carbon black or char- 
coal used was not specified. Consequently the excellent results were 
not duplicated here, although the following types of charcoal were used: 
bone, NucHar C-190N, and Darco G-60. In the Japanese article the 
whiskers were grown in an atmosphere of air plus the gases given off 
by the starting materials. 

The article concluded that the hydrogen reducing atmosphere was 
produced by reaction of carbon with water present in the charcoal 
followed by the growth reaction. 

H 2 0 + C — ► H 2 +CO 
liq FeX 2 + H 2 Fe + 2 HX 

The intention of the following series of experiments was to modify 
the reference experiment in order to explore the possiblity of using 
helium as a carrier gas for transportation of selected materials into 
the chamber in order to enhance whisker growth. 

1(a). Experimental 

3.0 g FeCl 2 • 4 H 2 0 + 1.0 g NucHar C-190 brand 
charcoal were ground together with a mortar and pestle and placed in 
a quartz boat. The boat was placed within a 16mm diameter quartz 
tube approximately 10 cm long between two smaller iron boats filled 
with Darco G-60 brand charcoal. Helium flow was 0. 15 CFH; time of 
reaction was ten minutes: temperature was 730°C. 
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He 



Fig. 3 Apparatus using carbon as a hydrogen source and 
heluim as a carrier gas. 



Results: A fair number of short (1-2 mm) whiskers were grown, how- 
ever, most of the metal was reduced as platelets. A white coating 
on many whiskers was also observed. 

1(b). Experimental 

The above experiment was modified by bubbling the He 
through two erlenmeyer flasks in order to moisten the atmosphere. 
Temperature = 730°C. Time = 10 minutes. He flow = 0. 2 CFH. 

Resul ts: The yield was very poor. No whiskers were formed; rather 
the metal was reduced ns platelets approximately 0.5 mm square. 



19 



It had been expected that the extra water in the atmosphere would 
increase the hydrogen concentration by: 

C + H 2 0 CO + h 2 
However, this was not effective. 

1(c). Experimental 

The above procedure was repeated with both charcoal boats 
upstream. He flow - 0. 15 CFH. Temperature = 730°C. Time = 30 
minutes. 

Results: A few whiskers were observed on the surface but none were 
found when the mass was probed. The appearance of the mass in the 
quartz boat ranged from powdery at the upstream end to granular at 
the downstream end. All whiskers grown appeared at the powdery 
upstream end. By placing both boats upstream of the starting material 
the hydrogen concentration was increased at that end as can be seen by 
the growth occuring at the upstream end. One may conclude the poor 
yield of experiment 1(b). was due to lack of hydrogen. 

1(d). Experimental 

3. 0 g FeCl 2 • 4 H 2 0 + 1.0 g NucHar C-190N charcoal was 
prepared as before and placed in a quartz boat. Directly upstream of 
the quartz boat was placed a small iron boat containing a ground 
charge of FeC^ • 4 ^0. Furthest upstream was another but longer 
iron boat containing a charge of Darco. Helium flow rate = 0. 2 CFH. 
Temperature = 800°C. Time = 11 minutes. 

Results: No whiskers were in either boat. A few platelets were 
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formed in the quartz boat. Most of the FeCl 2 • 4 H 2 O charge had 
vaporized, however, a few flakes of reduced iron were present. The 
expected results had been that the extra source of FeC^ vapor would 
increase the growth, yet this was not the case. 

Discussion: 

Using helium as a carrier gas was not advantageous. The 
helium flow, as was observed in experiment 1(a). , did not drastically 
reduce the whisker growth as compared to growth in air. This was 
predictable since the reactions occurred in the charcoal mass at vexy 
short distances between hydrogen sources and growth sites such that 
the low velocity helium didn't displace the hydrogen before the growth 
reaction occurred. 

There are two possible reasons for no enhancement of growth by 
introducting extra FeC^ vapor to the system as was done in experiment 
1(d). First, the hydrogen source was so small as to create a very 
weakly reductive atmosphere. The other possibility was that the FeCl 2 
must be in liquid from to react so as to form whiskers. Thus the 
FeCl 2 vapor did not react. 

An overall factor that would influence the form of the iron produced 
would be the super saturation of the halide vapor. It is generally be- 
lieved (Ref. 3) that whiskers grow in vapors whose supersaturation ratios 
are not far in excess of 1.0. Introduction of foreign material such 
as water vapor could affect that ratio such that the iron is produced in 
a different form. 
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2. Experiments with Ferric Oxide and Ammonium Bromide Mixtures 



Purpose: 

It had been previously shown in. this laboratory that a 
mixture of Fe 2 03 and NH 4 Br could produce excellent yields of iron 
whiskers under a hydrogen atmosphere. The reaction was thought to 
be: 

H 2 + 4NH 4 Br + Fe 2 C >3 4NH 3 + 2 FeBr 2 + 3 H 2 0 
FeBr 2 + H 2 -> Fe + 2 HBr 

An undesirable characteristic of this particular mixture was that 
the whiskers were fastened securely to the substrate. Consequently 
the crystals were not easily removed. The following series of 
experiments was executed to find a method of growing iron whiskers 
which could be detached from the substrate and be relatively free from 
flaws. Two methods were used to accomplish this end. The first 
involved a series of experiments which examined the addition of a 
loose substrate to the starting material. Carbon which reacted in the 
environment was added, and then various inert materials were added to 
accomplish this goal. 

The second method examined was the growth of whiskers on a 
surface not in contact with the starting materials. It was hoped that 
transporting the reactive vapors to sites not in contact with the original 
mass would lead to finding a site which would allow their removal. 

2(a). Experimental: 

To 11.0 g of ground NlI jBr, 7.0 g Fe 2 03 was added, and the 
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mixture ground together in a mortar. (Henceforth this mixture will be 
called 7-11.) 3. 0 g of this mixture was placed on an iron liner which 
was set into a quartz boat. The charge was reduced under a hydrogen 
atmosphere at a flow rate of 0. 5 CFH at 730°C. for 13 minutes. 

Results: A high yield of whiskers of varying sizes was obtained, the 
longest of which was over 1 cm. Growth often occurred in clumps 
along the sides of the liner. At the bottom of the liner the metal was 
reduced as a foil with some shorter whiskers growing out. Some of 
the whiskers were observed to have balls of material at their ends. 

This might have indicated a VLS mechanism. If this was the case then 
an impurity was present fulfilling the requirements of a nucleating 
site as stated in the previous discussion of VLS mechanism. 

2(b)l. Experimental: 

3. 0 g of 7-11 was ground with 1. 0 g of NucHar C-190N 
charcoal and the resultant mixture reduced with hydrogen. Temperature 
= 730°C. Hydrogen flow rate = 0. 5 CFH. Time = 13 minutes. 

Results: The mass was very powdery and the volume was greatly 
increased over that of the starting materials. A few whiskers were 
identifiable. When the mass was probed with a spatula, no whiskers 
were found below the surface. 

2(b)2. Experimental: 

The mass of charcoal was reduced to 0. 25 g and the 
experiment rerun under the same conditions. 

Resul ts: The mass was very voluminous with the surface having 
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a metallic nature rather than charcoal as was seen in experiment 2(b)l. 
Some short whiskers were entwined in the mass at the surface. No 
whiskers were found below the surface. 

2(b)3. Experimental: 

3 g of 7-11 was ground with 0. 1 g NucHar C-190N and 
reduced under the conditions of the previous experiment. 

Results: The mass more closely resembled that of experiment 2(b)l. 
(where no charcoal was added) than that of experiment 2(b)2. where 
0. 25 g charcoal was present. However, the yield was not as great as 
in experiment 2(b)l. Most of the whiskers grew on the liner sides while 
the metal was reduced as a foil at the bottom. Several grew in clumps 
and were often bent or misshapen. 

When these clumps were probed by a spatula carbon particles were 
found beneath the surface as well as a few whiskers. Beneath the foil 
at the bottom of the liner much carbon was found as well as a few 
metal particles. 

2(b)4. Experimental: 

Experiment 2(b)3. was performed using a charge of 3. 0 g 
of 7-11 and 0. 15 g of NucHar C-190N charcoal. 

Results: Two types of whiskers were formed. On the sides long, thin 
whiskers were present attached to the melt. At the bottom of the liner 
a mixture of charcoal and metal was present with short whiskers 
disbersed throughout. 

2(b)5. Experimental: 
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The experiment was repeated with 3. 0 g of 7-11 and 0. 12 g 



NucHar. 

Results: The boat was tilted during the reduction, consequently, there 

was a hard melt on the lower side. On the upper side some very long 
thin whiskers grew from a powder mass. Charcoal was found beneath 
this mass. 

2(c)l. Experimental: 

A 3.0 g charge of 7-11 mixture was placed in the liner, and 
0. 12 g of NucHar was spread over the surface. Otherwise the condi- 
tions were the same. 

Results: Three kinds of masses were present. The usual shiny 
metallic foil with carbon granules beneath had collected at the bottom. 
No whiskers were piesent here. Higher up the side of the liner a 
second type of mass was found having a gray powdery nature. This was 
a mixture of mostly carbon with metal. A few short whiskers were 
present. The third type of mass was situated above the second type on 
the side of the liner. A number of long whiskers were present and 
could be extracted relatively easily by using tweezers. A general 
observation was that the growth direction was perpendicular to the 
surface of the liner. 

2(c)2. Experimental: 

Experiment 2(c)l. was repeated. 3. 0 g of 7-11 mixture 
was placed in the steel liner, and a large amoung of NucHar was 
sprinkled on the surface such as to give a thin coating ever the entire 
surface. 
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Results: Overall, the entire mass was black with spots of metal 
coating on the surface. At the upstream end the mass was almost 
entirely carbon- like in appearance. As the mass went downstream it 
became more metallic in appearance. The yield of whiskers was not 
as great as that of experiment 2(c)l. with the number fewer, and the 
size shorter and wider. A few areas of powdery mass were present 
with some long whiskers growing out. When the entire mass was probed, 
much charcoal was found beneath along with some metal and no 
whiskers. 

2(d). Experimental: 

3. 0 g of 7-11 mixture ground together with 0. 25 g Darco 
G-60 brand charcoal. Temperature = 730°C. Hydrogen flow = 0. 5 
CFH. Time = 15 minutes. (Darco was a coarser more granular 
brand of charcoal. NucHar was a light, powder charcoal. ) 

Results: Most of the metal reduced was in the form of a foil in the 
bottom of the liner, while some clumps of metal were attached to the 
side. These small masses were hard and made up of carbon and metal 
granules. Many whiskers also grew from their surfaces, however, 
several whiskers were deformed. 

2(e). Experimental: 

3.0 g of 7-11 was ground together with 0. 12 g NucHar. 
The hydrogen was bubbled through water before flowing into the cham- 
ber. Time = 10 minutes. Temperature = 730°C. Hydrogen flow = 0.5 

CFH. 
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Results: At the upstream end there was a gray powdery mass as 
described in experiment 2(c)l. A few long whiskers as well as many 
short thick ones were present. At the downstream end was a metal 
foil where some long but deformed whiskers grew from carbon particles 
along the sides. 

2(f)l. Experimental: 

2. 0 g of 7-11 mixture, 1. 0 g FeC^ • 4 H 2 O, and 0. 12 g 
NucHar were ground together and reduced under hydrogen. Tempera-- 
ture = 730°C. Time = 12 minutes. Hydrogen flow rate = 0. 5 CFH. 
Results: The predominate mass was a shiny metal foil. Many long 
whiskers were growing from the sides where the substrate was more 
granular. Shorter whiskers were growing from the foil surface. Most 
of the carbon was found beneath the foil when the mass was probed. A 
few powdery, gray masses were present at the downstream end with 
whiskers growing from them. 

2(f)2. Experimental: 

Experiment 2(f)l. was repeated increasing the amount of 
NucHar to 0. 2 g. 

Results: The entire mass was a large powder "cloud" with a few short 
whiskers growing in the mass. 

2(f)3. Experimental: 

Experiment 2(f)2. was repeated several times varying the 
mass of NucHar only, with the following results: 

Results: a. 0. 15 g - A typical metal foil was formed. 
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b. 0. 17 g - One powdery mass with a few short whiskers 



was observed. 

c. 0. 16 g - Here was a powdery surface containing a few 
long whiskers. The longer ones grew along the sides of the liner, and 
beneath the powdery surface was a metal foil with charcoal underneath. 
Some of the whiskers were very thin and only noticed when a ball was 
seen at their tip. 

2(g). Experimental: 

A series of experiments was performed with 1. 5 g of 7-11, 
1.5 g FeCl 2 • 4 H 2 O, plus varying amounts of NucHar. Temperature = 
730°C. Time = 13 minutes. Hydrogen flow = 0.5 CFH. 

Results: 2(g)l. No charcoal added - a very large yield grew from 

the metal foil. In general the whiskers were very long and thick in 
diameter. The thickest growth (in number of whiskers per unit area) 
grew from the upper edge of the charge while the sparsest was along 
the bottom of the boat. Several clumps with whiskers growing out were 
on the liner sides. The whiskers were attached tightly to the substrate 
and not easily removed. 

2(g)2. 0. 2 g - a gray, powdery cloud with a few short 
whiskers was formed. 

2(g)3. 0. 12 g - a foil was along the bottom. Some long 
misshapen whiskers grew from fine granular masses on the sides. The 
masses were hard, and whiskers were not easily removed by tweezers. 

2(g)4. 0. 15 g - most of the metal produced was in the form 
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of a foil. The foil wasn't shiny; rather it was gray. Some very long 
whiskers grew from granular masses on the side of the liner, two of 
which were of such length that their growth was stopped when their tips 
touched the chamber wall. 

2(g)5. 0. 16 g - again most of the reduced metal was foil, 
however, present were a few fairly long whiskers which were easily 
removed from a gray mass situated along the upper edges of the boat. 

2(g)6. 0.16 - the charcoal was first heated alone for 10 
minutes under a hydrogen atmosphere at 730°C. Then the charcoal was 
mixed with the standard mixture of 7-11 and FeCl 2 • 4 H 2 O. The re- 
sultant product contained three types of masses. At the upstream end 
the mass was light gray and powdery with no whiskers. The bottom was 
covered by a metal foil while a darker metallic region was present at 
the edge of the foil as well as in the downstream end. Whiskers grew in 
this region but were not easily removed. 

2(h). Experimental: 

A series of experiments was run using 1.5 g FeCl 2 * 

4 H 2 O + 1. 5 g of 7-11 ground together. The variable in the series was 
the amount of flake graphite added to the above mixture. Tempera ture 
= 730°C. Hydrogen flow = 0. 5 CFH and time was 10 to 15 minutes. 
(Time was not an important parameter if all reductions were allowed at 
least 10 minutes to react. ) 

2(h)l. 0. 21 g - a very large yield of long whisker was 
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obtained. However, it was not as large a yield as that observed when 
no graphite was present. Most were grown on the sides of the liner. 

2(h)2. 0. 42 g - most of the metal produced was in flakes, that 
is, the metal had formed a coating about the graphite flakes. Many 
short whiskers were along the upper edge of the liner as well as a few 
long ones. They were fairly easily removed as compared to the pre- 
vious experiment. 

2(h)3. 0. 30 g - a larger number of whiskers were grown than 
in experiment 2(h)2. but less than 2(h)l. They were not as easily re- 
moved as in experiment 2(h)2. 

2(i). Experimental: 

A series of experiments was performed with the mixture 
of 1.5 g FeC^ • 4 H 2 O and 1.5 g of 7-11. This was mixed with 
several different inert materials in varying quantities in an attempt to 
find a substrate with which a good yield of whiskers could be obtained 
while allowing them to be easily removed. Temperature = 730°C. 
Hydrogen flow rate = 0. 5 CFH. Time = 12 to 15 minutes. 

Results: 2(i)l. 0. 16 g Nucllar + 0. 30 g quartz powder - one large 

gray, powder cloud was obtained, with a few whiskers along the upper 
edges. 

2(i)2. O'. 10 g of NucHar + 0. 30 g quartz powder - most of 
the metal was reduced as a foil, and some very long whiskers were 
growing in the direction of hydrogen flow. 

2(i)3. 0. 30 g quartz powder - the substrate was very 
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brittle and many whiskers were grown, however, these were securely 
fastened to the base. 

2(i)4. 0. 7 g iron filings - a good yield of long whiskers was 
present; the crystals were attached to a metal substrate and not easily 
removed. Many iron filings were found beneath the metal foil. 

2(i)5. 0. 7 g iron filings + 0. 16 g NucHar charcoal - one large 
black cloud with a few short whiskers was obtained, and most of the 
metal was reduced to a gray foil over the charcoal. When this foil was 
probed two layers were found beneath. The upper layer was gray in 
color with metal grains reduced in the mass. The lower layer contained 
less metal and was darker in appearance. The iron filings were not 
found at the bottom of the liner as was the case in experiment 2(i)4. 

Also, compared with experiment 2(g)5, where no iron filings were used, 
the effect due to charcoal was enhanced. In that experiment the mass 
was largely a foil, while here it was less metallic. 

2(i)6. 0. 7 g iron filings + 0. 10 g NucHar - most of the metal 
was in the form of a foil with a mixture of carbon and metal grains be- 
neath the foil. One cannot distinguish by eyesight if these metal grains 
were formed by chemical reactions or were the iron filings in the 
starting materials. A powdery mass was at the upstream end with 
many whiskers growing there. 

2(i)7. 1.0 g iron filings + 0. 12 g NucHar - a metal foil was 

formed with a light, gray, powdery mass on its surface with some 
whiskers growing out. Many of these whiskers were easily removed 
with tweezers. 
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Discussion: 



The previous sixteen procedures were performed with two 
purposes. First, it was to improve the yield of the 7-11. This was 
done by increasing the amount of reducible halide vapor with equal 
weights of FeCl 2 • 4 H 2 O and 7-11 mixture. The second purpose was to 
find a method of growing a good yield of whiskers on a substrate which 
would allow them to be easily removed. 

This purpose was not accomplished by use of inert materials 
added to the mixture. It was found that the addition of any of the inert 
substances tested decreased the yield of the original mixture . In 
general it was observed that these inert substances were denser than 
the FeC^ melt. Consequently they sank to the bottom of the liner and 
were often found beneath a metal foil which had accumulated there. 

Charcoal was an inert material as far as the growth reaction of 
iron whiskers was concerned. However, at 730°C it reacted by re- 
leasing water vapor. The amount of hydrogen produced by the reaction 
of carbon with water can be considered negligible compared with the 
hydrogen flow from outside the chamber. In general, there was a 
critical mass of charcoal that was added such that when less than this 
amount of carbon was added the result was a metal foil barren of 
easily removed whiskers. When an amount of charcoal in excess of the 
critical amount was added, the result was a voluminous mass with a 
few short whiskers. When an amount very near this critical value was 
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added as in experiment 2(g)5. , the whiskers were relatively easy to 
remove. However, die yield was very low. 

Taking an overall view of the entire series one can formulate 
several trends. The best whiskers grew on the liner sides. When the 
starting materials melted, growth sites were created as the edge of die 
melt retreated down the sides of the liner. The sites mentioned where 
whiskers were relatively easy to remove were due to carbon clinging to 
the sides and die growth occuring on their surfaces. 

3. Growdi of Iron Whiskers onto Various Surfaces. 

Since there was very limited success in growing easily remov- 
able whiskers by mixing various materials within the starting materials, 
another approach was followed. Here an attempt was made to separate 
the growth sites from the starting materials by transplanting the vapor. 

The first series of experiments dealt with the use of steel covers 
for the quartz boats. The reasoning here was that the vapor would re- 
duce on the covers, and the whiskers possibly be removed from there. 

By experimentation it was hoped that a substrate could be applied to 
these surfaces and would accomplish diis objective. 

The second series of experiments dealt in growing whiskers on 
various sharp metal projections. Various procedures were were used 
here to find how whiskers grew from tacks and pins in different posi- 
tions. 

3(a)l. Experimental: 

A charge of 7-11 mixture was placed in a steel liner which 
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was positioned in a quartz boat. The boat was covered with a flat 
steel plate with six, one- eighth inch holes drilled along the length of the 
plate. The cover was not air tight and would allow hydrogen to flow in 
under the edges of the plate as well as through the holes. Tempera- 
ture = 750°C. Hydrogen flow = 0. 5 CFH. Time = 30 minutes. 

Results: Most of the metal was reduced in the form of a foil on the 
bottom of the steel liner. A fair yield of whiskers was obtained, how- 
ever, it was inhibited as compared to a yield where the plate was not 
present. Many short whiskers were growing on the underside of the 
plate from very small sites which resembled crumpled bits of metal 
foil. High density of growth was found in and around the holes, while 
only one or two whiskers were found on the top surface around the holes. 
The whiskers attached to the plate were not easily removable with 
tweezers. 

3(a)2. Experimental: 

Experiment 3(a)l was repeated under the same condi- 
tions with two, one- eighth inch holes drilled through the plate. How- 
ever, the procedure was changed as follows: The chamber was purged 
with helium, and the boat placed in the hot part of the furnace for 15 
minutes. Then the hydrogen was allowed to flow at 0. 5 CFH. This 
was intended to have the starting materials melted and vapor starting 
to escape before the reducible atmosphere reached the boat. Then the 
vapor might be reduced while on the upper surface of the plate. 

Results: Results weren't drastically changed from experiment 3(a)! . 
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Fewer whiskers grew on the underside of the plate with the highest 
density growth around the holes. Whiskers were present in the boat 
although the yield was less due to the FeCl2 ’’boiling off” before the 
atmosphere became reductive. No whiskers grew on the upper surface 
of the plate. 

3(b)l. Experimental: 

2. 0 g of 7-11 and 1. 0 g FeCl 2 • 4 H 2 O were ground and 
placed in a metal liner. The metal liner was placed in a quartz tubing 
which was long enough to extend from the hot region of the furnace to 
the rubber stopper sealing the chamber off from the air as shown in the 
schematic. 




CONTAINING STARTING 
MATERIALS 

Fig. 4 Apparatus for experiments dealing in vapor transport. 



Helium was passed into the quartz boat as shown while hydrogen flowed 
in the atmosphere of the chamber. The boat was covered by a steel 
plate as shown with two, one- eighth inch holes drilled through it. I lhs 
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was done so that the vapor would flow through the holes onto the upper 
surface of the plate before coming in contact with the hydrogen. Temp- 
erature = 7 3(K>C. Time = 20 minutes. Hydrogen flow = 0. 5 CFH. 
Helium flow = 0. 2 CFH. 

Results: Most of the iron was reduced in the boat with a yield which 
was very good despite the conditions. A few whiskers had grown on the 
bottom of the cover but were difficult to remove. In the holes there 
was a high density of short whiskers. No whiskers were on the upper 
surface of the plate. A light colored film was present on this surface 
and stayed in the same condition over a period of days. Over the same 
period of time, however, the bottom of the plate corroded. 

3(b)2. Experimental: 

The above experiment was repeated with helium flow in- 
creased to 0. 3 CFH. 

Results: The results were the same as those of experiment 3(b)l. ex- 
cept the yield was less while the whiskers were of the same character. 
After three days the extreme ends of the plate had a corrosion product 
on the upper surface as well as over the entire bottom surface. On the 
upper surface, the middle section of the plate between the holes and a 
short distance toward the ends on either side was clear of the corrosion. 
The film reported in experiment 3(b)l. was in this area. 

3(b)3. Experimental: 

The experiment was repeated with the following changes: 

The steel liner was shoved as far up the tubular part of the quartz 
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boat as possible instead of directly beneath the holes as before. This 
was to prevent the hydrogen from reacting with the vapor at the liner’s 
surface so that it might excape to the upper surface of the plate. 

Helium flow = 0. 5 CFH. Hydrogen flow = 0. 4 CFH. Time = 30 minutes. 
Results: No whiskers grew in the area covered by the tube, although 
a few did grow at the downstream end which was beneath the plate. 
Growth occurred on the upstream end of the bottom surface of the plate. 
Whiskers also grew in the upstream hole but not in the downstream one. 
A whitish film was present on the upper surface but no whiskers. 

3(b)4. Experimental: 

Experiment 3(b)3. was repeated with helium flow at 0. 6 
CFH and hydrogen flow of 0. 2 CFH. 

Results: No whiskers were in the steel liner under the tubular part. 
Very few were at the downstream end of the liner as well the underside 
of the plate. A few were in the upstream hole but none in the down- 
stream hole. 

3(b)5. Experimental: 

Experiment 3(b)4 was repeated with the only change being a 
small amount of Fe 20 ^ was sprinkled on the upperside of the steel 
plate to provide possible whisker growth sites. 

Results: Results were the same as those of the previous experiment 
with no whiskers grown on the top of the plate. 

3(c). Experimental: 

The conditions of experiment 3(b)5 didn't provide growth 
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sites on the upper surface for any possible escaping FeCl 2 vapor to re- 
act to form whiskers. Since growth sites occur on the bottom of the 
plate the following procedure was carried out: a liner containing just 

a charge of FeCl 2 • 4 H 2 O was set into the quartz boat. A steel plate 
was placed upside down over the liner as shown. 




Fig. 5 Alteration of apparatus for experiment 3(c) 



The quartz boat was placed in the hot part of the furnace for ten 
minutes under just a helium atmosphere. It was expected that the 
starting material would' "boil" and possibly deposit some nucleating 
sites upon the plate. 

Results: The expected small pieces of foil commonly seen on die under- 
side of the plate in previous experiments were not present. Rather, the 
area was discolored possibly due to a film of material FeCH. One can 
say the small pieces of foil were deposited only when the materials 
were under a hydrogen atmosphere where the growth reaction was 
occurring. 
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3(d)l. Experimental: 

Two tacks were positioned in a steel liner with their tips 
sticking up out of the charge of 3. 0 g of FeC^ * 4 H 2 O. The liner was 
placed in an open quartz boat and reacted under just a helium atmo- 
sphere flowing at 0. 5 CFH. Time = 10 minutes. 

Results: Many whiskers were grown on the liner walls along a line 
defined by the original limits of the charge's volume. A very dense 
growth was on the rodlike part of the tack. The whiskers grew in a 
radial direction like the spokes of a wheel growing from the hub. There 
was a sharp line of demarcation on the rod where whisker growth 
stopped, approximately where the surface of the original charge touched 
the vertical part of the tack. 

3(e). Experimental: 

A tack was placed in a charge of FeCl 2 0 4 H 2 O as before 
in experiment 3(d)2. The boat was placed in the hot part of the furnace 
(730°C) under a helium atmosphere. >-Flow was 0.5 CFH. After 10 
minutes the boat was withdrawn, and a film on the sides of the liner 
was observed as well as on the tack. Small crystals (not whiskers) were 
also on the body of the tack. A new tack was situated in the liner next 
to the other one. The boat was replaced in die furnace for 10 minutes 
with a hydrogen flow of 0. 5 CFH. 

Results: The whisker yield was poor but in die same areas as that of 
the previous experiment. Whiskers appeared on both tacks, but growth 
was more extensive, and the crystals longer on die tack that had been 
under both helium and hydrogen atmospheres. 
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Discussion: 



Acceptable growth on the upper surface of the plate was not 
obtained. There were indications that FeC^ did, however, reach that 
surface and was reduced there. One may be able to conclude this from 
the noncorrosive film appearing near the holes. Iron whiskers in the 
boat and underside of the plate did not corrode with time. However, the 
steel liner and plate, which were not pure iron as were the whiskers, 
corroded. Thus, one may conclude the film on the upper surface was 
pure iron produced by reduction of FeCl 2 . If this was true, then the 
situation was one where the reduction occurred in whisker form on one 
side of the plate and as a film on the other side. 

Several reasons could be responsible for this phenomenon. The 
supersaturation ratio of FeCl2 above the plate could be a value not 
conducive to whisker growth, or growth sites for whiskers may not be 
present on the upper surface. On the underside, whiskers grew from 
small pieces of material present due to the reaction of the starting 
materials. 

This might be a future technique for studying nucleating sites. A 
series of experiments would have to be conducted such that starting 
materials would be subject to reductive conditions for a small interval, 
removed, observed, and replaced. A progression of nucleating site, 
whisker, and termination should be observed. This observation should 
be easier than trying to observe growth in the boat where sites are 
often covered by a film of metal. 
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The tack experiments seem to confirm that whisker sites tend to 
be above or follow the edge of the region bounded by the starting mater- 
ial. Whisker growth tended to be on the sides of the liner and on the 
upper part of the tack. This could show that a liquid film in the growth 
area cannot exceed a certain thickness, or it will poison the growth site. 

Experiment 3(e) showed that the nucleating sites are not formed by 
just the initial heating of the starting materials, rather hydrogen may 
be necessary. Sites for growth could be found even after 10 minutes of 
heating FeCl2 at 730°C under a helium atmosphere. 

B. GROWTH OF COPPER WHISKERS 

In reference 7 a method of growing large numbers of copper whis- 
kers was discovered. The materials in the experiment were CuO, CuCl, 
and graphite mixed in' weight ratios of 100:10:5. Experiments were 
carried out in this article using a furnace similar to that used here to 
grow iron whiskers. The above starting materials were placed in a 
boat which was set inside a 16mm diameter quartz tube between two 
pieces of graphite. This tube was then fitted inside the larger diameter 
tube making up the furnace chamber. Both ends of the chamber were 
sealed by inserting solid rubber stoppers. Thus, the growth of whiskers 
occurred in an atmosphere of air and those gases accumulated from the 
reaction, CO and CO2 in particular. 

Pure copper in whisker form was thus produced by reduction of CuO. 

2 CuO + CO - ^ C112O + CO2 
Cu 2 0 + CO-» 2Cu +CO2 
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CuCl was also a source of copper metal reduced by hydrogen pro- 
duced from water in the graphite. 

H 2 0 + C -> CO + h 2 
2 CuCl+ H 2 — >2 Cu+ 2 HC1 
2 HC1 + Cu 2 0-» 2 CuCl + H 2 0 

This method of whisker growth was of special importance due to 
the reported very high yields. Copper was said to have been produced 
in whisker form at nearly 100 percent of the theoretical yield. 

1(a). Experimental: 

3. 0 g CuO (powder), 0. 3 g flake graphite, and 0. 15 g CuCl 
were ground together and reduced as described above. These amounts 
were the weight ratios of 100:10:5 as prescribed by the reference article. 
The apparatus was modified as follows: (1) Instead of graphite, iron 
boats were filled with Darco brand charcoal and used to supply the re- 
ductive atmosphere. (2) The gas was allowed to escape into a water 
trap upon expanison. Time = 20 minutes. 

Results: No whiskers were detected. The surface of the mass was 
black and the interior contained red particles which were probably un- 
reduced Cu 2 0. On the metal liner was a copper colored film covering 
almost the entire surface. 

1(b). Experimental: 

6. 0 g CuO (wire form), 0. 6 g flake graphite and 0. 3 g CuCl 
were ground together and used as starting materials. Temperature = 
730oC. Time = 20 minutes. 
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Results: The results were the same as those of experiment 1(a). Much 
of the CuO hadn’t been reduced to CU 2 O which was present as red grains. 
One could determine this because CuO was still in the distinguishable 
pellet form in which it had been before heating. 

1(c). Experimental: 

3. 0 g CuO (powder), 0. 3 g NucHar C-190 brand charcoal, 
and 0. 15 g CuCl were the starting materials. Temperature = 730°C. 
Time = 20 minutes. 

Results: The mass resembled one homogeneous, voluminous ’’cloud" 
throughout. Looking with the naked eye the appearance was one of 
copper powder. When seen under the stereomicroscope at 30 power 
one was able to see a network of very small whiskers entwined through- 
out the mass. The mass was held together such that when probed with 
a spatula, the entire "cloud" would move. 

1(d). Experimental: 

3. 0 g CuO (powder), 0. 3 g Darco, and 0. 15 g CuCl were 
the starting materials. Darco was a coarser, more granular charcoal 
as opposed to the powdery NucHar used in 1(c). 

Results: The mass closely resembled that of experiment 1(c). and that 
reported in the literature (Ref. 7). Some "balls" containing many 
whiskers were observed. When penetrated, the mass was seen to 
assume more of a charcoal nature with depth. 
l(d)2. Experimental: 

3. 0 g CuO (wire), 0. 3 g Darco, and 0. 15 g CuCl were 
reduced under the previous conditions. 
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Results: The cloud was more voluminous and of a more copper color. 



No charcoal was on the surface as was observed previously. Whiskers 
were below the surface, however, as the mass was probed deeper the 
yield of whiskers decreased and the mass assumed more of a char- 
coal character. 

l(d)3. Experimental: 

3. 0 g CuO (powder), 0. 48 g NucHar, and 0. 15 g CuCl were 
reduced for 20 minutes at 730°C. 

Results: A yield no higher than that of l(d)2. was obtained. The mass 
was of a harder, cakey nature than before. 
l(e)l. Experimental: 

3. 0 g CuO (powder), 0. 3 g Darco, and 0. 15 g CuCl were 
reduced at 700°C for 41 hours. 

Results: There was high whisker yield, but the mass was essentually 
of the same voluminous nature. 

1(e) 2. Experimental: 

Experiment l(d)4 was repeated at 600°C. 

Results: The yield was the best obtained thus far, and it was easier to 
detect the strands of whiskers by the naked eye. The whiskers had a 
shinier copper color compared to those grown at temperatures in ex- 
cess of 700°C. Still, the whiskers are entwined and difficult to separ- 
ate. 

l(e)3. Experimental: 

Experiment l(e)l was repeated at 500°C for 22. 5 hours. 
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Results: No perceptible improvement or. decrease in yield was obtained. 



l(e)4. Experimental: 

Experiment l(e)l was repeated at 400°C for 12 hours. 
Results: Very poor whisker growth was observed, but, some very 
small diameter whiskers were observed such that they were only noticed 
when light refelected from them. Only a small percentage of starting 
material reacted. When the surface was disturbed by probing more 
whiskers could be detected by the reflection of light. 
l(e)5. Experimental: 

The experiment was repeated at 336°C for 48 hours. 

Results: A low yield was obtained, however, some very small diameter 
whiskers which were easily removed were observed. These were 
longer than those obtained at higher temperatures. 
l(e)6. Experimental: 

The temperature was raised to 430°C and the same mixture 
of new material was reduced for 18 hours. 

Results: A typical entwined, high yield was obtained. Many "balls" of 
tangled whiskers were observed. A few whiskers were long, straight 
and untangled such that they could be removed by the use of tweezers. 
l(f)l. Experimental: 

3. 0 g CuO, 0. 6 g Darco, and 0. 15 g CuCl were placed in a 
furnace at 750°C under standard procedures for one hour. However, 
the usual two accompanying boats of Darco were not used to produce the 
reductive atmosphere. 
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Results: Many isolated areas of tangled whisker growth were present 
on a surface predominately charcoal in nature. A few longer, straight 
whiskers were present. When the surface was disturbed more whiskers 
were observed when their movement reflected light. The growth below 
the mass was of the same nature. 

1(h). Experimental: 

A charge of 3. 0 g CuO, 0. 3 g Darco, and 0. 15 g CuCl was 
prepared and reduced under the following schedule: 
l(h)l. 15 hours at 330°C. 
l(h)2. 11 hours at 370°C. 

1(h) 3. 14 hours at 540°C. 

After each period of reduction, the inner quartz tube containing the 
quartz boat and two iron boats was removed, and the growth observed. 
New Darco was placed in the iron boats between each period, but the 
same charge of starting materials was used throughout. This was done 
to find how the growth progressed as the mixture was heated. 

Results: l(h)l. A low yield of tangled whiskers was observed at 

each end of the boat. A few short, straight whiskers were in the mid- 
dle section of the mass. Possibly the tangled growth could have been 
due to the proximity of the charcoal boat to these areas. 

l(h)2. The yield became greater. When a draft of air was 
directed at the surface, one could see small diameter whiskers by light 
reflection. Otherwise, they were difficult to see. These whiskers were 
seen in the middle section which was beginning to acquire areas of 
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tangled growth. Both ends of the mass were covered by whiskers en- 
twined in clumps. The small diameter, straight whiskers could have 
been present but obscured by the other type of growth. 

1(h) 3. The mass was of the typical high yield appearance 
expected at this temperature. Clumps of tangled whiskers were ob- 
served over the entire surface, while the small, strai ght whiskers 
could not be observed. 

l(i). Experimental: 

l(i)l. 3. 0 g CuO, 0. 3 g Darco, and 0. 15 g CuCl were re- 
duced for 10 hours with no end boats containing charcoal. The temper- 
ature was the lowest used in the growth of copper whisker at 316°C. 

l(i)2. The very same mixture was heated at 380°C for 21 
hours after its appearance had been noted. 

l(i)3. The appearance was noted, boat placed in the furnace 
again, and reduced at 730°C for 2. 5 hours. No end boats were used for 
extra charcoal. 

Results: l(i)l. A few short but straight whiskers were grown. Two 

relatively long ones were present on the sides of the liner. Whiskers 
growing on the surface were easily removed by tweezers. 

l(i)2. There was little more whisker growth. Some tangled 
clumps appeared in the middle section. 

l(i)3. The tangled growth increased, but the yield was still 
low. Some red crystals were present indicating Cu 9 0 that had yet to 
be reduced to Cu. No increase in short, straight whiskers was 



observed. 
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l(j). Experimental: 

A fresh mixture of the materials used in experiment l(i) 
was reduced without end boats at 730°C for one hour. 

Results: The typical high yield of entwined copper whiskers was ob- 
tained. 

Discussion: 

The reaction did not occur using graphite as reported in the 
literature, however, both types of charcoal used were successful. 
Perhaps the Japanese authors made a mistake in translation. Best re- 
sults were obtained at a temperature of 600OC. 

The high yields reported in the literature were duplicated. How- 
ever, these yields were of a tangled nature, not easily separated, or 
of the desired straight and thin appearance, it was possible to obtain 
this growth at lower temperatures, but at the same time, yield was 
drastically reduced and time of reaction increased. 

The concentration of the reductive gases, H 2 and CO, was of im- 
portance. At lower temperatures the production of these gases would 
decrease while at the same time the straight thin whiskers appeared. 
Also, at lower temperatures tangled growth appeared at the ends of the 
mass which would be in closer proximity than the center where the de- 
sirable forms appear. This could allow one to predict that in a highly 
reductive atmosphere tangled, unrestrained growth occurs, while the 
more perfect whiskers appeared in less reductive atmospheres. This 
is a common generality of crystallization where the most perfect 
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crystals grow at slower rates. As added evidence for this idea, the 
electron scanning miscroscope showed very few of the perfect whiskers 
in mixtures grown with iron boats containing extra charcoal. If a mix- 
ture of tangled growth and perfect whiskers was present and undetected 
under the visual microscope, the high magnification of the instrument 
should have detected many of these perfect whiskers. 

A suggested series of experiments to detect the role of CuCl in the 
mechanism will now be proposed. The effect of the amount of CuCl 
added should first be determined by varying the amount present. Since 
the reduction of CuCl depends upon the production of hydrogen obtained 
from the reaction of water and carbon, one would suspect that the small 
amount of water present would quickly react or evaporate away at 800°C. 
Thus the hydrogen supply is limited and appears only in the initial 
stages of the experiment. Therefore, could this mean that the reduction 
of CuCl is the initiating reaction producing the growth sites? To prove 
this, it might be productive to react the starting materials under a 
hydrogen atmosphere with carbon that had been heated previously under 
helium to remove water. The only reaction occurring should by the 
production of Cu by reducing CuCl. The mass should be examined 
carefully to see if this is so. Then reduce this same mass under an 
atmosphere of air and try to detect if growth was initiated from the 
sites produced by CuCl reduction. 
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IV. COLD TIP THEORY 



A. POSTULATED MECHANISM 

The VLS mechanism has had sound experimental evidence to prove 
its validity. However, this mechanism was applied only to nonmetals 
due to the lack of a suitable low melting alloy to serve as the condensing 
site for the iron whiskers. 

The screw dislocation theory has never been conclusively proven to 
be the mechanism for growing pure metallic whiskers. Much contro- 
versy has surrounded the second mechanism. A screw dislocation 
would be very difficult to detect due to its very size, also, this dislo- 
cation would be very difficult to detect under the theory that growth 
stops when the dislocation disappears. Field emmission microscopy 
(Ref. 8) and electron scanning microscopes have to date failed to show 
conclusive evidence of screw dislocations. 

The screw dislocation mechanism for metallic whiskers has had 
extensive mathematical treatment by Burton, Cabrera, and Frank in 
Ref. 4. Mathematics of diffusion controlled whisker grown from the 
vapor also has been given considerable attention (Ref. 9, 10), while 
heat flow in whiskers has been relatively ignored. The cold tip mechan- 
ism, as postulated in this paper, is valid only for endothermic reactions 
occurring on the whisker surface. As the reaction proceeds, heat is 
drawn from the material. The point of highest curavature, which is at 
the tip, will have the most heat withdrawn per unit volume. 
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Consequently, the temperature of the tip will be the coolest part of the 
whisker. 

The reaction material is adsorbed on the surface of the whisker 
or a liquid film creeps up from the base and diffuses to the region of 
lowest energy on the whisker. In the screw dislocation mechanism, 
this region of lowest energy is the dislocation. However, a cold tip 
would serve the same purpose by removing energy f rom the diffusing 
material and thus decreasing its ability to translate to a different part 
of the whisker or evaporate to the atmosphere. There is a higher pro- 
bility of finding reaction material at the tip, thus the reaction is more 
likely to occur in that region. Since the reaction is endothermic, more 
heat will be withdrawn from the tip and the tip will be cooled to an even 
greater extent. This phenomenon of a colder and colder tip is referred 
to as thermal runaway. 

A mechanism where the heat flow is the important consideration 
has been proposed previously to explain dendritic growth from solutions 
(Ref. 15). Due to the heat of crystallization, there is a heat flux into 
a cubic crystal. The corners and to a lesser extent the edges are able 
to dissipate the heat more easily than the flat surfaces. Consequendy 
the sides of the cube will accumulate heat faster than the edges and 
corners. This will be an important factor if the material has a low 
thermal conductivity. Otherwise this phenomenon of geometry will be 
offset by the heat being rapidly distributed throughout by conduction. 

For the case of low conductivity, growth will be stimulated at the 
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comers and to a lesser extent the edges. Thus, growth will be rapid at 
the eight corners and approach zero at the faces of the crystal. This 
complements the cold tip mechanism very well where the growth is even 
more directionalized due to the growth reaction being endothermic 
rather than exothermic as in the dendritic growth. 

Since the reaction must be endothermic, one is able to analyze the 
nature of the reacting materials at the growth site. For iron whiskers 
growing at a typical growth temperature of lOOOoK the reaction is: 

FeCl 2 + H 2 (g)-» Fe(s) + 2 HCl(g) 

Where the phase of FeCl 2 is unspecified to make the following 
points: (Ref. 16) 

(1) If FeCl 2 is in the gaseous phase, the reaction is exo- 
thermic by approximately 6,000 calories per mole. 

(2) If FeCl 2 is liquid, the reaction is endothermic by ap- 
proximately 24,000 calories per mole. 

(3) If the above reaction were to occur on another area of 
the whisker, and an iron atom formed in the liquid phase, it may then 
migrate to the tip. The reaction at the growth site would then be exo- 
thermic due to the heat of crystallization of Fe. 

All three of the above cases could occur within the framework of 
screw dislocation theory, but case two is the only one which would lead 
to a cold tip. In addition, the first case where the FeCl 9 vapor reacts 
directly at the tip is rendered less likely by the equilibrium constant; 
at 1000°K the reaction is to the left. A similar reaction occurs 
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for copper whiskers and the three cases correspond to those shown 
above. 

With FeCl2 in the liquid phase, other effects of the cold tip can be 
noted. If the FeCl2 tends to condense on the vertical part of the whis- 
ker before diffusing to the tip this region is heated by the heat of vapo- 
rization along the sides thus decreasing the chance or reaction in this 
area. The tip is thus colder relative to the sides. A colder tip will 
cool the atmosphere around it condensing FeC^ which releases the 
heat of vaporization destroying or at least decreasing the temperature 
of the tip. The tip would become less attractive as a reaction site 
compared to the sides. Consequently the whiskers would react unifor- 
mly over the surface until the cold tip was restarted by the endother- 
mic reaction and the heat flow due to the high curvature at the tip. 

A more probable prediction based on the FeCl2 condensation 
would be that an equilibruim state would be reached such that the heat 
flux into the tip by condensation and convection; equal heat removed by 
chemical reaction. The temperature at thermal equilibruim would 
still be substantially less at the tip compared to the sides of the whis- 
ker to insure preferential growth at the tip. 

A different scheme of transport of liquid FeCl2 postulates that the 
material diffuses from the base where it is present as a melt. The 
liquid is "hot" relative to the metal in that it is nearly the temperature 
of the substrate. As it advances, the sides of the whisker are heated, 
and the tip becomes cooler relative to the sides. This would tend to be 
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dominating effect when the whisker is short. As the whisker becomes 
longer the liquid from the base has further to travel. The whisker near 
the substrate would approach the temperature of the base due to its dis- 
tance from the cold tip and proximity to the base. Consequently the 
liquid film would tend to evaporate from the sides once the tip has 
grown a certain distance from the base unless the atmosphere was 
supersaturated. This would be one explanation of why the atmosphere 
must be at least slightly supersaturated for whisker growth. 

Summarizing, the effects of transport of liquid gives the following 
insight into whisker growth by the hypothesized cold tip mechanism. 
Growth is initiated by a nucleating site on the substrate. This site 
does not necessarily have to be a cool area but must be an area condu- 
cive to crystallization. Initial growth is due to'liquid from the melt 
diffusing up the whisker sides. An intermediate period of growth occurs 
with a cool tip at a thermal equilibrium between the endothermic react- 
ion and exothermic condensation. Finally, growth is terminated when 
the atmosphere is less than saturated and vaporization of the liquid 
occurs faster than it is supplied. This terminating situation would 
occur when the whisker has grown above the halide vapor or when that 
vapor has simply been depleted from the starting materials. 

A very rough mathematical approximation for the temperature 
gradient can add a degree of rigor to the qualitative predictions. Sup- 
pose iron was formed by the liquid FeC^ reacting with hydrogen over an 
area of one square centimeter and grew a crystal from that surface at 
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100 JK / sec. In 100 seconds, 1 cubic cm would have been deposited at 
a heat of reaction of approximately -3500 calories. If this heat was 
supplied through just the tip by conduction from the metal one could 
predict: 

q = -k A 1 X 

H * x 

Ceil 

q = rate of heat transfer. (- 3500 XUO sec) 

Ceil 

k = thermal conductivity. (0. 16 Q ) 

S6C \_>i cm 

\ r r Op 

= temperature gradient. 
o x cm 

A = area (1 cm^) 

NT 

Solving for : 

^ T = 219 °C 

ax rnf 

Qr, over 100^/f the temperature difference between base and tip is 
2. 2 °C. This figure is very approximate and represents the greatest tem- 
erature gradient possible taking into account just heat transfer entirely 
through the tip. However, this calculation shows the tip could be sign- 
ificantly cooler than the rest of the whisker. 

B. THE MATHEMATICS OF HEAT FLOW 

1. Differential Equations to be Considered 

Before the whisker can be modelled, it is necessary to exam- 
ine the differential equations describing heat flow. A very clear ex- 
planation of basic heat flow equations is given in reference 11. 

For pure conduction in any material, the heat flows with the tem- 
perature gradient -jX- from an area of higher temperature to a lower 



56 



temperature region This is expressed by Fourieir's law: 

u } T 
^ = - kA Tx~ 

Where: q = Heat transfer rate. 

k = Thermal conductivity. 

A = Area. 

T = Temperature, 
x = Distance coordinate. 

In three dimensions the equation for the unsteady state condition, 

that is the change of tempera ture with time is: 

* 2 T + a 2 T + + q - pc 

ix 2 by 2 ^ Z 2 k k 77 
Where: q = Energy generated per unit volume. 

= Density of the material.- 

c = Heat capacity of the material. 

Since the problem under consideration does not have a heat generat- 
ing factor the equation becomes: 

v 2 t _ 1 ^T 
^ at 

Where: c< =~ ^T6 ~ and is called the thermal diffusivity 

for a specific material. 

Another important relationship used in heat transfer is Newton's 
law of cooling due to convection: 

q = hA(T w - r k> ) ’ 

Where: h = Film conductance whose units are heat per 

unit time, area, and degree of temperature. 
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T w = The temperature of the material at the 
boundary. 

Too = The temperature of the fluid flowing by 

the solid boundary at a distance such that 
its temperature is unaffected by the 
solid's temperature. 

h is not a constant for any material. It depends on geometry, 
velocity, and other variables and must be determined experimentally 
for each circumstance. 

Another type of heat transfer needed to be examined for the pro- 
blem under consideration is radiation. The net transfer of energy be- 
tween two black bodies at temperatures Ti and T 2 is described by the 
Stefan - Boltzmann law; 

q = cr A (Tj 4 - T 2 4 ) 

Where : <T = The Stefan - Boltzmann constant. 

This law is only valid for perfect emitters, therefore, to describe 
"gray" bodies compensating proportionality constants are used. 

2. Finite Difference Techniques 

In the analysis of heat flow, the temperature of a point depends 
upon all the points in its universe. Several complicating factors are 
found in the growth of whiskers in relation to the heat flow. The geo- 
metry of the boundary, as well as the boundary conditions, and the 
movement of the boundary due to deposition of material make the sol- 
ution very difficult by analytical techniques. Consequently, the problem 
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is more conveniently solved using finite difference techniques (Ref. 12). 
This method substitutes discrete values to represent the continuous 
variable in its assigned region. In the problem under consideration, 
it is desired to relate a temperature with its position on the whisker. 

To do this a grid was superimposed over the shape and temperatures 
assigned to the grid points. The convention for the grid on the whisker 
was as follows: the-central axis of the rod was called the y axis. The 
whisker was assumed to have cylindrical symmetry with the x axis 
representing the radial distance from the center. Thus the plane exam- 
ined in this paper was that defined by the central axis and a line perpen- 
dicular to that axis extending to the surface of the whisker (Fig. 6). 

With a grid now superimposed, one may derive suitable relationships 
called difference equations for the mathematical analysis. 

Calculations involving finite differences depend on replacing the 
differential equation with a difference equation. A difference equation 
approximates the differential equation by relating a set of points on the 
grid in a manner specified by the differential equation. This set of 
points is called a molecule and the differential equation to be solved 
in this paper is Fourier's law of heat conduction. 

Some basic finite difference equations which are relative to the 
differential equation can be seen from Figure 7 where A is the unit grid 
size. 
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Fig. 8 

Whisker and its grid overlay 
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Fig. 9 Finite difference relationships 



These equations evaluated at T 0 are: 

^T = T 3 -Tj 
Tx 2A 

a 2 T = T 3 - 2T +Tj 



d T = T, - T 2 



2A 



d T = T 4 - 2T n +T 



0 



C) X^ A -1 A 

However, these relationships only hold for grids of constant A . 
In the shape under study, a complicating factor arises when the irreg- 
ular boundary intersects the lines of the grid between points. A differ- 
ence expression for the first and second partials of temperature with 
respect to nonuniform distances can be derived as follows: 

Consider the grid with the boundary as shown. 

Where: A = the size of the unit grid. 

PjA= the distance from point 0 to point j. In the 
above case it can be seen P| and P 9 are 1.0 
P 3 and P^ are fractions. 
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Fig. 10 An irregular boundary and grid. 



The differential equation contains the expressions 



YT 



^ 2 T 

c) y2 



A 

and — — — for which substituting expressions are derived as follows: 

^ z 2 

First expand T(x 0 +A, y 0 ) and T(x 0 , y 0 +B), about (x 0 , y 0 ) 



which are the coordinates of point 0, by use of a Taylor series. 
T(x 0 + A, y G ) = T(x 0 , y G ) + A-^- T(x 0 , y Q ) 

+ . . . 



+ ~ 2 ~ TT2 T(x °’ y ° ] 



B" 



<} 2 



T(xq, y G +B) T(xo, Yq) +B 3 y T(x 0 , ^ ^( x o?yo) 

+ - - . 
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Substituting A =£3 A and then A = -P^A into equation one as well 
as B = P4A and then B = - P2A into equation two, the following 
equations are obtained after truncating the rest of the terms in the 



series. 



T( x 0 + ?3 A , y 0 ) = T(x 0 , y Q ) + P 3 A -jz T(x q , y Q ) + 

T (x 0 > y Q ) 

T(Xo -Pl A , y Q ) = T(xq, y Q ) - Pi A -|y T(Xo,y 0 ) + 



(P4A y 






T(xq, y Q ) 



\ 2 

T(Xo, y 0 + P 4 4 ) = T(x 0> y 0 ) + P 4 A T(x 0 , y G ) fy 

t(xq, y 0 ) 

^ (P4 A ) 2 

T(Xq, y Q - P 2 A ) = T(x 0 , y Q ) ~^T T ( x o> y 0 ) + 2 5^2 

T(xq, y Q ) 



Solving the first two equations simultaneously and setting Tj = 
T(xj, yj) the expression for the second partial of temperature in the x 
direction at point 0 is derived. 

<^ 2 T 0 _ 2 1 



x 2 A 2 R? + p l P q L 



_ 2 



y t q 

^ y 2 A 2 P|+P 9 P 



2 4 



To - T 0 (l+Zl)+ II T, 



T 4 -T 0 (n|±-)+ii T 2 



( 1 ) 



( 2 ) 



b 2 T 



An expression for — can be found by examining Laplace's 
equation in terms of r, e, y where : (Ref. 13 ) 
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CL> 



x = r cos © 



z = r sin © 



y = y 



<) 2 t ^ 2 t ^ 2 t 

^ x 2 ^ y 2 c) z 2 



A 2 t 1 At 1 

" r t r 1 9 

A r 2 r ^ r 



A 2 t [ ^ 2 t 

b ©^ 3 y2 



y 




Due to the cylindrical symmetry of the whisker r = x = z, and the 

problem can be analyzed in just the x - y plane where © = 0. Laplace's 

equation thus reduces to : 

A 2 T _L 

c) x c) x 

Substituting this relationship into Fourier's law, the conduction 
equation becomes: 
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( 3 ) 



b 2 T <) 2 T 1 ^ T = 1 T 

<}x2 <)y2 x^x 0(^t 

Difference equations for the first partial of temperature with re- 
spect to x can be derived by simulanteously solving the same Taylor 
series expansions used to derive equations (1) and (2). This results in: 



1 JT 

X ^ X 



p l P 3 

x A (Pi + P 3 ) 





(4) 



Notice that for points not adjacent to the boundary the fractions P^, 
?2> P 3 and P 4 are 1.0. Equations (1) and ( 2 ) reduce to: 

^ 2 T _ T1 - 2T n +T3 <i 2 T _ T? - 2T n + T4 

ix 2 = A 2 by 2 " A 2 

and equation (4) reduces to: 

^ T T 3 - T x • 

Xx 2 A 



For the study of unsteady state heat the following convention in 
nomenclature is followed: 

Tj = Tempera ture at point j. 

Tj, + 1 = Temperature at point j at one time interval 

later. 

Using Fourier’s law of conduction assuming no heat source within 
the element: 





£ T 

yr 
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the temperature one interval can be calculated. 



AT = , - 2 ^ . 

7 rr v T 

AT = (X ( V 2 T) At 

Substituting the finite difference equation for V where the point 
under evaluation is not adjacent to the boundary: 



AT = T 0 , + 1 -T 0 = c* At 



Tl - 2T 0 + T 3 T 2 - 2T 0 +T 4 1 



T 3 - T 1 



X 



(5) 



T, 



o, +1 



= c< A t 



t 1 + t 2+ t 3+ t 4 T 3 -T x 
+ 



x A 



4<x A t 

+ To (1 - 7-J)(6) 



The term 



c< A t 

time increment such' that: 

A 2 



imposes a restriction upon the grid size and 



cX At 



= M > 4 



This factor is called the modulus, M, and if its value is less than 4 
this leads to an unstable condition. A nonrigorous illustration will show 
why this is so. Suppose that Tj, T2, T3, and T4 are equal and cooler 
than T 0 . Equation (6) reduces to: 

4Ti cx At 4 cx A t 

T^ . + T 0 (l- -— ) 



0, +1 



4 <x At 

A ~1 (Tl " T o> + T o 



M 



(Tl - T 0 )+T ( 



o 
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If: 



M > 4 




M= 4 



M <4 



The condition MO predicts that Tp has cooled below its surround- 
ings, or, in other words, energy has flowed "uphill". One would imtui- 
tively guess that T Q would approach at steady state. When M = 4 
the system is conditionally stable, and when 4 the system is stable. 
Thus, once a A is chosen the range of At is fixed. The most accur- 
ate results for finite difference techniques is obtained for small A 's. 
However, A t must be very small to maintain stability, and the steady 
state is approached slowly. If one desires to obtain the solution more 
quickly, he must increase At, and, at the same time A , thus losing 
grid accuracy. If a boundary is superimposed on the grid such that for 
points adjacent to that boundary, Pg and are less than 1.0, the mod- 
ulus is different than than calculated for the points on the interior. Here 
A is smaller, thus, ©< At must be decreased to maintain stability increas- 
ing computation time on the computer to reach steady state. For in- 
stance, in the whisker analysis, it was necessary to adjust of At such 
that M = 4 at the points adjacent to the edge, but for interior points the 
sameckAt gave M = 20 since the effective grid size changed. 



The equation: 



Ti ,i = T. +<X( V 2 T) 4 t 

h ' 1 J 



Where: V^T = the sum of equations (1), (2), and (4), can 
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calculate the new temperature for points interior to the boundary pro- 
vided proper values of Pi, ? 2 , P3 and P4 are used. However, for 
points on the boundary the image method is convenient to use (Ref. 14). 

The image technique involves overlaying an imaginary grid line 
perpendicular to the point under consideration. 




Fig. 12 

Image point technique. 

One is able to simulate the desired flux across the boundary by 
working with image temperatures Ii and I2. Ii is the actual tempera- 
ture at a distance one A inside the whisker normal to the point under 
consideration T2. The value of 1^ can be found by interpolating between 
temperatures which are on grid points. For the above situation these 
temperatures would be T 00 , Tq^, T ,Tq _ p T_jq, and T^, and the 
method used for interpolation will be shown later. 

12 is an imaginary temperature whose value is manipulated to 
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obtain the desired flux. For example., an insulated boundary is one 

through which heat does not flow. Therefore, the temperature 
^ T 

gradient is 0 where x is understood to be the direction normal to 

0 " 

the surface. In difference equation form this is: 



! 2 - ; i n 

Jx = 26 



Therefore: 

i 2 = I i 

To calculate T 2 + ]_, equation (6) must be used. This assumes the 
boundary is linear over the interval from point 1 to point 3. The four 
temperatures T]_, Tg, ^ and and their corresponding values of P^, 
^2’ ^3» anc * ^4 wou lb be used. 

Interpolation to find I 2 involves many steps, the first of which is 

determining the slope at point 2. The method of calculation assumes a 

quadratic function describing the positions of points 1, 2, and 3, where 

the coordinates of these three points are: 

(Xx.Y^ 

Fig. 13 

Points used in approximation of slope 




The function can be linearized about (X 2 , yg) by Taylor series. 



, by 

?1 ~Y2 + cix 



1 b^y 



2 (Xl “ X 2 } + 2 dx2 



(X]_ - x 0 ) 2 + . . . 
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2 



, dy 
73 = ^2 + dx 



i d 2 y 

(x 3' x 2> + 5^2 



(x 3 - x 2 r + 



Truncating after the third term and solving for ^ evaluated at 
point 2: 



dy = 1 

3^ 2 • x 3 - x 2 - (x 3 - x 2 ) 2 

< X 1 - x 2 ) 



y3 ‘ y 2 + 



( 7 ) 



It is necessary to use the above slope to calculate the physical pos- 
ition of in relation to Tqq in order to interpolate a value for 1^. 



\ 




Fig. 14 

Determining position of image points. 



From Fig. 14- one can derive the distances XT and YT which are 
used later in the interpolation process. 

Ax^ + 2 



- 1 2 

ay= Ax(-j±r 
dx 

dx 
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Substituting: 



Ax 2 + flx 2 ( 1 ) 2 = A 2 

dy_ 
dx 

= A'< 1 + 1 1 > ^ (8) 

(dy } 2 

dx 



By a similar derivation: 



y 



i + d^ 2 
( dx> 



) 2 



Then: 

xT = X2 - Ax - 
yT = y 2 - * y - y x 

1-^ was determined by using a six-point interpolation formula in the 
following manner: (utilizing the points shown in Fig. 15) 




Fig. 15 

Points used for six point interpolation. 
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A second order, two-dimensional Taylor series yields: 



h = T 00 + < T x>x = 0 < xT > + i < T xx>x = 0 < xT > 2 + <Vx = 0 <y T > (9) 
y = 0 y = 0 y = 0 



+ i <T yy ) x = o (y T > 2 + i ( T xy>x = o < x r> (yT) 

y = 0 y = 0 

Where: 

= 0 “ \ y eva l uate( l at Tqq. 
y = 0 

^ T 

(T y ) x = 0 = Ty“ evaluated at T 00 . 
y = 0 

, and so forth. 

Referring to the basic finite difference expressions for first and 
second partials one can see them in terms of the unit grid size: 



< T x>x = 0 

y = o 
(T xA = o 

y = 0 

( T y) x = 0 
y = 0 



T m- T-10 

2 



- T 10 - 2T 00 +T. 10 



T 01 ' T 0-l 



(Tyy) x = 0 T 0 1 - 2 Tqq + Tq_ 

y = 0 

An expression for (T X y) x _ q was formulated by solving expan- 

y = 0 

sions where temperatures are known: 

Tqi = Tq + (Ty) x = o + k (^yy)x = 0 

y = 0 y = 0 

Tfi - Tq + (T x ) x _ q + 2 (T xx ) x = 0 (^y)x = 0 2(Tyy) x = o 

y - 0 y = 0 y = 0 y = 0 

+ k (Txy) x = 0 
y = o 
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Solving simultaneously: 

(Txy)x = 0 = ^ (T^i " % - (T x ) x = q “ 2 (T xx ) x = q ) 
y = 0 y = 0 y = 0 

Substituting T x and T X y.: 

( T xy)x = 0 “ ^11 ~ 2 T 0 i - 2T^q + 2Tqq 

y = 0 

Thus, Ii may now be calculated by substituting the proper values 
of xT, yT and the above difference equations into equation (9). 



C. COMPUTED RESPONSES OF THE WHISKER TO BOUNDARY 
CONDITIONS 



When a whisker is in the growth environment there are several 
boundary conditions to be met at the whisker surface. From the inter- 
ior heat is transferred to the surface by conduction. 



q = -kA 



^ T 



x 



The heat flow leaving the surface is governed by several effects. 
Convection and radiation depend upon the temperature difference be- 
tween the surface and atmosphere. Both factors are complicated by 
the film of adsorbed gas or liquid present on the boundary as well as 
the flow of hydrogen past the surface. Thus, the proportionality con- 
stants for each effect would be very difficult to determine experiment- 
ally. 



However, a few predictions may be made from some simple calcu- 
lations. It has previously been shown (pg. 5"(?) that the flow of heat for 
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an iron whisker growing by reduction of FeC^ at 100^/{/sec is approx- 
imately -35 cal/sec due to the endothermic reaction. Since convection 
as well as radiation depend upon temperature differentials, an initial 
approximation was made where the chemical reaction was the only fac- 
tor taken into consideration. 

An additional factor that is present in the cold tip theory of whisker 
growth is the heat imparted by the diffusion up the sides of the whisker 
by the liquid halide. However, this was not treated here due to the 
complexity of this particular type of heat transfer. 

Two types of boundary conditions were imparted to the surface and 
the responses studied. First, as a control to see the type of heat flow 
within the geometry of a whisker, the surface was insulated. The calcu- 
lations used to study this involved giving the surface a thermal shock at 
time = 0. Then the whisker responded by a heat flow due only to con- 
duction from the base. This response is shown by the isotherms of the 
whisker drawn on the next few pages. 

A word of explanation concerning the isotherms is appropriate at this 
time . The figures are scaled to contour the temperature arrays from a 
computer, consequently, the horizontal scale is different from that of the 
vertical. Temperatures of the isotherms are on the right side with the 
value in the upper right hand corner indicating the temperature at the tipi 
An isotherm of constant value is drawn across the base to indicate the 
the heat source. The isotherms resembling steps are due to the fact that 
only three significant figures were stored in the arrays. For example, a 
temperature of 29.95°C is stored as 30.00QC. 
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Radius (microns) 



Fig. 16 

Initial condition of pulsed whisker 
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Radius (microns) 



Fig. 17 

Isotherms after 100 iterations. 
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22.2 
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30.0 



- Fig. 18 

Isotherms after 200 iterations. 
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Radius (microns) 

' Fig. 19 

Isotherms after 300 iterations. 
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20.9 
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25,3 
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Fig. 20 

Isotherms after 400 iterations. 



79 




Radius (microns) 



Fig. 21 



Isotherms after 500 iterations. 
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The time shown was calculated from the interval c< a t. To under- 



stand this refer to equation (5): here o( A t was a quantity used to mult- 
iply the finite difference. In the computer analysis, a value for ex' A t 
of 0.003 was determined to be the largest number allowing a stable 
solution, that is, a modulus approximately equal to 4. o<for iron at 
1000°K is 0.04-^HL. or 4 x lO^^Li^/sec. Since <*At equals 0.003 the . 
true time interval per iteration is 0. 75 x 10“ ^ seconds. 

The thermal shock was applied at the surface in the following 
manner: a temperature difference of 30°C was assumed to take place 
between the base and surface. The temperatures at all other points 
were calculated as a function of their distance from the surface. This 
function was: 

T = 30 • (1 - exp (-B* D 2 )) 

Where: D = Distance from the surface. 

B = An arbitrary number describing the attenuation 
of the shock with distance into the interior. 

From an examination of the positions of the isotherms with time 
one can see that the tip is the coolest point on the whisker. It can be 
seen after 500 iterations the temperature is dependent almost exclusi- 
vely upon the distance from the base as is expected. 

The next series of isotherms was drawn from the treatment of the 
chemical reaction at the surface. It was assumed that the whisker and 
atmosphere were at the same temperature and the surface was covered 
uniformly by the necessary halide for the growth mechanism. At 
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Fig. 22 

Whisker isotherms after 100 iterations. 
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Whisker isotherms after 200 iterations. 
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Whisker isotherms after 300 iterations. 
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Fig. 25 

Whisker isotherms after 400 iterations 
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Fig. 26 



Whisker isotherms after 500 iterations. 
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time = 0 the reaction was started. The reaction was treated as if it 
occured uniformly over the surface. Also, it was assumed that it 
withdrew enough heat at all points to cool the surface by 0. 1°C per 
iteration. The surface was again assumed to be insulated. 

From an examination of the series, one can see the tip is still the 
coldest point on the whisker. In comparison to the previous case where 
the endothermic reaction was not under consideration, one can see the 
entire whisker is gradually cooling instead of warming. 

D. SUGGESTIONS FOR FURTHER ANALYSIS 

The computer modelling done to date on this problem has been 
limited in its consideration of FeCl 2 oh the side of the whisker. One 
can see in the series of isotherms calculated with consideration of the 
chemical reaction (Figs. 22- 26) temperatures are low adjacent to the 
sides of the whisker. This situation would be eliminated by taking 
into account heat transfer between the diffusing liquid and metal. At 
the base the liquid is at the temperature of the base. As it moves the 
temperatures inside the metal are lower than the base due to the cold 
tip. Consequently there is a temperature gradient present such that 
heat flows into the metal. Temperature of the metal increases and the 
liquid decreases. Thus, as the tip is approached the liquid has been 
cooled appreciably. 

The second consideration is the kinetics of the rate of diffusion 
with distance from the base. Growth rate will decrease due to the 
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greater distance and hence longer time necessary to migrate to the tip. 
This is an area of study deserving extensive consideration and needs to 
be applied to the present model. 

V. SUMMARY 

The work of this paper can be summarized as follows: the cold tip 
mechanism put forth could be the correct mechanism for the growth of 
metallic whiskers. At the least, an examination in this direction has 
shed new ideas on how to explain whisker growth phenomenon. It pre- 
dicts the diffusion of the metal halide in the liquid state from thermo- 
dynamic considerations. What is needed to continue the theory is an 
extensive analysis of the diffusion of this liquid keeping the cold tip in 
consideration. 
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APPENDIX A 



Computer Modeling 

Due to the particular boundary under consideration, an extensive 
explanation of how the mathematics of section B were adapted to the 
computer is necessary. First it was necessary to make several basic 
assumptions about the system. The whisker was assumed to have the 
dimensions of 10 microns in length (L) at the stage of growth under 
observation. The physical boundary was described arbitrarily by the 
equation: 

Y = L - (AX 2 ) / (r 2 -X 2 ) 

Where: 

• • Y = The vertical coordinate coincident with the 
longitudinal axis. 

X = The radial distance from the center. 

A = A dimensionless parameter which changes 
the desired curvature of the whisker. The 
larger the value of A the sharper the curva- 
ture. 

r = The radius of the whisker at full growth. 

Since it was decided that the whisker was perfectly round, the 
above equation described a parabolic cone. This symmetry simplified 
the problem, and thus it was only necessary to treat the points on one 
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side of the central axis. A unit grid size of 0. 25 microns was laid 
over the whisker in order to apply finite difference techniques (Fig. 6). 
Several assumptions were then made for the environment of this whis- 
ker. The base was assumed to be resting upon a heat source, thus the 
base was kept at constant temperature. No attempt was made to account 
for the cooling of the atmosphere by the endothermic reaction, that is, 
it was considered insulated from the whisker surface. Thus, the 
temperature of the atmosphere was constant everywhere it contacted the 
metal surface. Also, the whisker boundary was not allowed to move as . 
if it were growing. This greatly simplified the problem and confined 
the study of heat flow to just the curvature of the boundary and the cool- 
ing due to the chemical reaction. 

First it was necessary to formulate a two dimensional array of 
temperatures called T(I, J). I corresponded to the element in the x 
dimension and J in the y dimension. The grid size of 0. 25 microns was 
called A in the derivations for the difference equations. Here it was 
given the name DEL. Thus, the physical position (X, Y) of the (1, J) 
element which was interior to the surface element was: 

x = (I - 1) • DEL 
y = (J -1) • DEL 

Complications arose when it was desired to represent the surface 
points. The convention used can ixtoreieasily be understood from the 
following diagram: 
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Fig. 27 

Assignment of temperatures to grid points 



The above diagram shows the two situations that arose. Tempera- 
tures were assigned to points where the boundary intersected the grid. 
If the boundary intersected the horizontal grid line of value J then this 
temperature was stored in the (1+ 1, J) element. If the boundary 
intersected the vertical grid line I shown, its value was stored in the 
(I + 1, J + 1) element. 

To determine how the boundary intersected the grid a test for each 
point and its relation to the boundary was formulated. The vertical 
distance from a point to the surface was called YPY. 

YPY = YP - Y 



91 



Where: 

YP = L - (AX 2 ) / (r 2 - X 2 ) 

Y = (J - 1) • DEL 

The horizontal distance to the surface was XPX. 



XPX = XP - X 



XP = 



( L- Y) r 2 
. (A + L- Y) 



1 

2 



X = (I - !)• DEL 




Fig. 28 

Testing element (I, J)’s proximity to surface 



If XPX was less than DEL the boundary intersected the horizontal 
grid line between the I and I + 1 element. Consequendy the surface 
temperature at this intersection was stored in the (I + 1, J) element. 
Similarly, if YPY was less than DEL, the temperature of die point 
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represented at the intersection of the vertical grid line and boundary 
was stored in the (I + 1, J + 1) element. 

If one recalls the difference equations (1), ( 2 ), and (4), it can be 
seen that for all points not adjacent to the surface the fractions P 4 , P2, 
P3, and P^ were such that they were all equal to 1.0. However, for 
the adjacent points^Pg and P 4 became fractions while P^ and P2 remained 
1.0. When a point was tested in the program and found that either or 
both XPX and YPY were less than DEL, P3 and P^ were calculated and 
stored in the one dimensional arrays P3O) and P 4 (I) respectively by: 

P 3 (J) = XPX / DEL 
P 4 (I) = YPY / DEL 

These arrays were one dimensional since PgQ) and P 4 (I) were 
unique for each horizonal and vertical grid line as seen in Fig. 8. 

At the same time, it was desirable for later logic statements to 
assign a number to the various points in the array. A new two 
dimensional array IN was formulated such that: 

IN(I,J) = 1 for all points interior to the surface. 

IN(I, J) = 2 for all points on the surface. 

IN(I, J) = 0 for all points outside the whisker. 

The computer logic necessary for this testing and assigning 
technique can be seen under subroutine ASSIGN in the appendix. 

The next phase in adapting the whisker to the computer model was 
to assign the values T 4 , T2, T3, and T 4 to evaluate T Q ^ + 4 by equation 
(6). The general case .for points inside the boundary was: 
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Where: 



ZO = T(I, J) 

Z1 = T(I-1, J) 

Z2 = T(I, J-l) 

Z3 = T(I + 1, J) 

Z4 = T(I, J+l) 

ZO = TO, Z1 = Tl; and so fourth. 

However, exceptions to the above case had to be made for several 
cases. If 1=1, that is, if the x coordinate = 0, Zl= T(I+1, J) due to the 
symmetry of the whisker. Along the base^ which was a heat source at a 
constant temperature called TINT, Z2 = TINT.. For points directly 
below points where the boundary crossed the vertical grid lines 
Z4 = T(I + 1, J+ l)'since the value of the temperature at the intersection 
was stored in the (1+ 1, J +1) element. The computer logic for these 
steps is contained in subroutine TNEW. 

Once the above values were assigned the change in temperature per 
unit time was calculated using equation (6). This change for every point 
interior to the surface was assigned to the two-dimensional array 
DELTAT. The logic for this operation is contained in subroutine 
DERIV. 

Now that the change in temperature per unit time for interior 
points had been calculated it was necessary to treat the edge points. 
First it was necessary to assign N - Y coordinates to these surface 
points. 
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Fig. 29 

Assignment of coordinates to surface points. 



In the above diagram, one can see the usefulness of the IN array 
whose values are . shown at the grid points. When a point was recog- 
nized to have a value of 2 its x coordinate was assigned: 

x= (1-1 +P 3 (J)). DEL 
y = (J-i) • DEL 

These assignments were valid for points like A and C. For points 
of type B: 

x = (1-1) • DEL 
y = (J-2+ P 4 (1-1) ) • DEL 

The next phase of the problem was to calculate the value of the 
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image point one DEL normal to the surface at the point under evaluation 
First the slope of the function at that point was determined. This was 
done by evaluating equation (7). However, several different express- 
ions were needed to find the three sets of coordinates (xj., y^), (X 2 , y 2 ) 



and 


(x 3 , y 3 ). (See Fig. 30) 






The expression for these coordinates were: 


(a) 


x 2 = (1-2+ P 3 (J) ) . DEL 


y 2 = 0-1) • DHL 




X 1 = (1-2 + P 3 (J + 1) ) • DEL 


yj = ] • DEL 




x 3 =(l-2+P 3 (J- 1))' DEL 


y 3 = 0-2) • DEL 


(b) 


Same as (a) except: 






x x = (1-2) • DEL 


yi = 0-1 + ) • del 


(c) 


x 2 = (1-2) • DEL 


y 2 = 0-2 + P 4 (I-l)) • DEL 




xi = (l-3+P 3 (J-l)). DEL 


yj = 0-1)- DEL 




x 3 = (1-2+ P 3 (J-1) ) • DEL 


y 3 = 0 - 2 )- del 


(d) 


Same as (a) except: 






x 3 = (1-1)* DEL 


y 3 = 0-2+ P 4 (l) ) • DEL 


(e) 


x 2 = (l-2+P 3 (J) ) • DEL 


y 2 = 0 - 1 ) ; del 




x x = (1-2) • DEL 


y 1 = (J-l + P 4 <l-l))- DEL 




x 3 = (1-1) • DEL 


y 3 = (]- 2 + P 4 (l) ) • del 


(f) 


Same as (c) except: 






x x = (1-3) • DEL 


y 1 = (J-2 + P 4 (I-2))- DEL 




The final values left to calculate to use equation ( 6 ) were the 



proper values for the fractions of DEL: Pj_, P 2 > P 3 , and P 4 • Due to 
the image points being at a distance of one DEL, P^ and P 3 equaled 
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Different cases for assignment of 
coordinates used for approximation 
of slope. 



Fig. 30 
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1.0. ?2 a nd were calculated by the distance formula between any 
two points in a plane. 

P 2 = ( (x 2 - X3) 2 + (y 2 - y3> 2 )i / DEL 
P 4 = ( (xj - x 2 ) 2 + (yi - y 2 ) 2 )i / DEL 
Using the appropriately calculated values equation (6) was applied 
to the edge, and the change in tempreature over one time interval was 
calculated for those points. 
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oonononooooonnoooonnnononooooooonooooononoonoooooooooonoonooonnoonnnoono 



THIS PROGRAM IS DESIGNED FOR THE STUDY OF HEAT FLOW IN 
A WHISKER OP A SHAPE DESCRIBED BY THE FUNCTION; 

Y = E L L- ( A* X** 2 ) / ( R**2-X**2 ) 

WHERE • • • 

A = A * PARAMETER WHOSE VALUE IS PROPORTIONAL 
TO THE DEGREE OF CURVATURE 
ELL=LENGTH OF THE WHISKER 
R=RADI US OF THE WHISKER 



•*r *r- 'V' n** -t- • a r “T- *r- *r- t- ^ nr -v nr «r* -t' ■n- nr -t- ■n' -nr *r* *tt ^ -Tr* •nr •nr *v» *t~ -t- "p *v *¥• *7* t* ■tt ~r- t* *r* -t* *tt *tt ■'?> *■»- ■t* -t 



DATA CARD ONE. CONTAINS THE FOLLOWING INFORMATION TO BE 
ENTERED IN A FORMAT (3 I 5 , 11 F10 . 5t 1 1 5 ) IN THE FOLLOWING 
ORDER: 

N=THE NUMBER OF GRID ELEMENTS IN THE X DIRECTION. 
M=THE NUMBER OF GRID ELEMENTS IN THE Y DIRECTION. 
INTER8=THE NUMBER OF ITERATIONS DESIRED BETWEEN 
PRINTOUTS OF THE TEMPERATURE ARRAY. 

ELL=LENG TH OF THE WHISKER. 

TSURF= IN I TI AL SURFACE TEMPERATURE. 

TINT=TEMPERATURE OF THE BASE AND ATMOSPHERF. 
DEL=SIZE OF THE UNIT GRID. 

F=THERMAL CONDUCTIVITY OF THE WHISKER MATERIAL 
DIVIDED BY ITS HEAT CAPACITY 
RHO=DENSI TY OF THE WHISKER MATERIAL. 

T I M I NT =T I M E REPRESENTED BY ONE ITERATION. 

SI GMA=A PARAMETER USED TO INDICATE THE FLUX 
THROUGH THE BOUNDARY. 

F I NCR= A CONSTANT REPRESENTING THE TEMPERATURE 
CHANGE DUE TO CHEMICAL REACTION PER 
.ITERATION 

I DCI D=A NUMBER USED TO EXERCISE THE CONTOUR OPTION 
1=USE OF OPTION TO CONTOUR TEMPERATURES. 
0=DES I RE NOT TO CONTOUR. 



U# sJU JU JU w ^ sb ^ J/ *JU sL- «JU •JU *JU JU ^ JU JU *JU JU Jx -JU st> JU JU «X» «JU JU OU JU *JU ^ *JU JL* *JU *JU JU *JU J- ^ 

nr -y* '¥• -¥• t t v ^ nr n- *r ^ ^ nr nr nr nr nr nr n^ nr nr nr nr n^ nr n' nr nr nr nr n' nr nr nr n^ — p nr n* 



DATA CARD TWO CONTAINS A NUMBER CALLED TAU IN A 
FORMAT ( 1F10. 5) . THIS VALUE IS USED TO LABEL THE 
PRINTOUT INDICATING THE NUMBER OF ITERATIONS 
PERFORMED UP TO THAT TIME. THE PROGRAM WILL TERMINATE 
WHEN IT NO LONGER HAS A TAU TO READ. 






<JU <JU J* 4- 4# 4» JL> J. 4* ">*' 4' 4» 4< 4L> J-- 4> 4» 4» 4. 4^ <X <X> <JU 4-> V. 4 . 4* iV 4U <X> JU >JU 4. JU «JU 4> 4f 4^ JU JU 

<v ■t- “T* *r* 'r* *r* n- '/* •r* "1” - 4 * -r- *r* -t* *r -t* «v -v nr nr -r* *r t- nr ~p *1' ^ *r- o'* t* ^ -r* -n -c “t* -c 



DATA CARDS THREE AND FOUR CONTAIN THE TITLE FOR THE 
GRAPHS IF THE ARRAY IS CONTOURED. DO NOT INSERT THESE 
CARDS IF THE CONTOUR OPTION IS NOT EXERCISED. 



JUJUO^ 4 . 4 . JU JUJU ^ JU4«S<^ JU 4 . J.JU J. J. 4.JUXJU 4. J<>ki, 4 f 4U JUJU JU JU JUJU 4. .u. JU 4 . .U4. X J. X 4. 4> J.JUJ. JU 4 . 4. 4<4r' x ‘ 

T T T T' *1' «Y* <T* ^ *r- •T*' *t'“ •’T - . T* "T *r» -I"* «■,% ^ -I' -t- -I' •?* •'5 V *T' "T* ^ "i'' 



THE PROGRAM USED TO CONTOUR THE TEMPERATURE ARRAY IS 
CALLED CONTUR AND IS CONTAINED IN THE REFERENCE 
LIBRARY OF THE NAVAL POSTGRADUATE COMPUTER CENTER. 



■* Of OU JU 0< O# U# 4# 4^ JU 4* «JU 4- 4- hU JU 4> «L «JU» 4U JU *JU 4^ 4U *JU JU JU «JU JU JU *JU J- wU *JU 

r n^ n? nr* nr nr nr nr n*- nr nr nr n*- nr n*- nr nr nr n % ^r n^ nr ^r nr o' # r n' nr nr nr nr **p nr nr 
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c 

c 

c 

c 



c 



c 

c 

c 

c 

c 

c 

c 



c 

c 



c 

c 

c 

c 

c 

c 

c 



THE FOLLOWING ARRAYS ARE DIMENSIONED SUCH THAT FOR 
AN ELEMENT (I,J), THE MAXIMUM I IS AT LEAST N+2 AND 
THE MAXIMUM J IS AT LEAST M+2 

REAL*4 T(20,50) ,P3(50) ,P4( 20 ) , DE LT AT ( 20 , 50 ) ,TXX( 50,20) 
INTEGERS I N ( 2 0 , 5 0 ) 

THESE ARRAYS DEAL WITH CONTUR 
LOGICAL*]. LTG ( 3 ) 

RE AL*4 CL ( 6 ) 

REAL*8 T I TLE ( 12 ) 

EQU I VALENCE (TXX, DELT AT ) 

1 READ( 5,1 00) N ,M, INTERS , ELL, R, TSURF, TINT , DEL, F, A, RHO, 

*T I M I NT , S IGMA , F INCR , I DC I D 
100 FORMAT (315 ,11F5.2, 115) 

105 READ (5,110) TAU 
110 FORMAT ( IF 10 • 5 ) 

N1=N+1 
M 1=M+ 1 
M2=M1+1 
WRITE (6,115) 

115 FORMAT ( 1H1, T4, 'N', T9, *M* ,T13 , • I NTER8 • , T25 , • ELL • ,T35, 

* ' R ' , T 40 , ' TSURF ' ,T52, 'TINT* ,T61 , ' TI MI NT • , T73 , 'DEL' ,T79, 
*• SIGMA' ,T90, ' F I NCR ' , T100, • TAU' ,/ ) 

117 WRITE(6,120)N,M» INTER 8 ,ELL , R , TSURF , T I NT ,TI MI NT , DEL , 

* S I GMA , F INCR, TAU 
120 FORMAT (IX, 315,1 IF 10. 5,/) 

^ ’JUh <JL« %L» Of OL- sLh ^JU J- ^ O# O# JU Of OU JL> O- OU sb «Jb> O- sl» OU 0> vL OU Of Of <JU OU OU %Jb- OU O. O- Op O* O* 

^ - 7 ** ^ ^ -p 7^. 7^ /j-* ^ O* P* T* -H* * 7 * 'T* v 7** Jp #p ^p Pp ^ p* p 'p ^p *"p Jp 7^ -*p >p 7% >p- 



ASSIGN INITIAL TEMPERATURES 

OLh O' O- OU O# O# O' o^ O* OLh O# Oi* O^ Of vb JU sb sb O^ JU vb O' O' O' sb Ob O# O' O# O/ sb sb Ob O' O^ *Jb» Ob Ob sb «JL~ Or «b Ob O' Ob Ob O' O' Of 
/p pp -p ^ -p «p ^p -7^ rp ^p >p -tp p #p #p p, p. p o* t* p *r p p p p p p p p p p p p p p p p p p p p p p p p p 



CALL ASS IGN ( N, M, T, ELL, R, TSURF, TINT, DEL »XPX , YPY , P4 » P3 , 
* I N , A ) 

THE FOLLOWING DO LOOP WRITES THE ARRAY AS PICTURED 
WITH THE BASE AT THE BOTTOM OF THE PAGE 
NC 0UNT =-l 

121 NC0UNT=NCGUNT+ 1 
J=Ml-NC0UNr 

WRITE (6,130) (T( I, J) , 1=1, N1 ) 

130 FORMAT ( '0S2X18F6.1) 

IF(MCOUNT.EQ.M) GO TO 122 

123 GO TO 121 

122 WRITE (6,135) 

135 FORMAT (1H1,T4, 'THE IN(I,J) ARRAY IS',/) 

NCCUNT =-l 

124 NC0UNT = NC0UNT+1 
J=M1-NC0UNT 

WRITE (6,140) (IN(I ,J) ,1=1 ,N1) 

140 FORMAT (/13Io) 

IF ( NCGUNT . EQ.M ) GO TO 202 
GO TO 124 



=£ £ * * # sjs s$: if £ * 5{c * # sjt ajs sjc * :js s}: # * * * £ $ * £ * s{c ^ # 3jt # # # 



»*< »■> J* ■JL* sL> <JL< JU JU \U 

*1' -O' •O' O' O' '1' O' O' O' 'O 



ITERATE THE T(I,J) ARRAY THE SPECIFIED NUMBER OF TIMES 






S** 'U «J. J* v? •A' 4 jf 4 sU sU JU vL> vU JU J. JU X 

O' o' or or or O' '1' »r ■V' -r o* ^ *r- 'r *r* or or o'* •Y' o' -r* 



# # ?*C # * # # 5*C # # £ # 



-JL. .Ju. -X- -JU. X- — -J- 
o' o' -t- or or or or 



700 DO 136 I = 1 , I NT ER8 

90 0 CALL TNEW ( N, M »DEL, TI NT ,TSURF ,T , IN, DELT AT , P4 , P3 ,T I MI NT, 
*F , C , RHO , ELL, R, SIGMA, FI NCR) 

136 CONTINUE 

WRIT E ( 6 , 1 15 ) 

WRITE(6, 120)N,M, INTER8 , ELL , R , T SURF , T I NT , T I M I NT ,DE L , 
*SIGMA,FINCR,TAU 
NC0UNT=- 1 

126 NCOUNT =NC0UNT+1 
J=M1— NCOUNT 

WRITE (6, 130) (T( I,J) ,I=1,N1) 

IF ( NCOUNT .EQ.M) GO TO 125 



100 



ooooooo ooooooo no o on on nonoDon 



127 

125 

202 



600 

50 



141 



142 
190 

143 

144 

145 



950 



10 



146 



150 



100 

110 



GO TO 126 

READ (5,110) TAU 

I F ( IDCID.EQ.0) GO TO 700 

###*# 3j! A >jc 3* Jjt £3!; 

CONTOUR THE ARRAY OF TEMPERATURES 

READ 15»50) TITLE 
FORMAT (6A8) 

THIS DO LOOP EXTENDS THE SURFACE TEMPERATURES 
OUT INTO THE ATMOSPHERE FOR A NEATER GRAPH 
DO 145 J = 1 » Ml 
DO 144 1=1, N1 

IF ( INI I , J) .EQ. 1) GO TO 144 
IFUN1 It J) .NE.2) GO TO 143 
T SAVE = T ( I t J> 

GO TO 144 
T ( 1 1 J I =T SAVE 
CONTINUE 
CONTINUE 

THIS DO LOOP ROTATES THE T < I , J > ARRAY ON ITS SIDE 
MAKING THE RESULTING PLOT EASIER TO UNDERSTAND 
DO 10 1=1, N1 
DO 10 J=1 , Ml 
L=M2- J 

TXXIL, I ) =T ( I , J ) 

THESE VALUES ARE USED BY CONTUR 
LTG( 1) =.TRUE. 

LTG(2)=. FALSE. 

LTG(3 )=. FALSE. 

MX=5 0 
NL=— 6 
I W=5 
IH=7 

INSERT THE CALL CARD AS SHOWN IN THIS SAMPLE 
CALL CONTUR(TXX,Ml,NlyMX,CL,NL,TITLE, IW, IH,LTG) 

GO TO 700 
END 



SUBROUTINE ASS I GN( N, M , T , ELL , R , TSURF, T I NT , DEL ,XPX,YPY, 
#P4 , P3 , IN, A) 

REAL*4 T (20,50) ,P3(50) ,P4(20) ,DELTAT(20, 50) 

INTEGER*2 IN(20,50) 

RS=R*R 
M1=M+1 
N1=N+ 1 
M2=M+2 
N2=N+2 



t** : 



j «JU U» vU U# • 



U U< <JU 



-r* n' <r* *v* 



<JL> JU <JU X X X X X X X X X X X X JU 
T T T> O' T ^ T 'f* T"T"r *T* W V 



INITIALIZE ALL T(I,J) AND IN C I , J ) 



X X X X- X X X X 

n 5 -r* -r *r -r -v* 



X X X X X X X# X X 
— «' I ~ ' — 



uxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 
•* *¥* 'v 'T* -t* *v -t 5 * ■" -v -v “V ■*,* *v -¥* -t* **6 -r- ■T' -V 1 -n* «t' -t» *ir "T* 'T* <*v* 'n 



DO 110 J = 1 , M 1 
DO 100 1=1, N1 
I N ( I , J ) =0 
T ( I , J ) =T SURF 
CONT INUE 
CONTINUE 



TEST PROXIMITY TO SURFACE AND ASSIGN P3,P4,IN 






101 



OOOOOOOOO 



DO 1000 J = 1 1 mi 
Y= ( J— 1 ) ^DEL 

XP=SQRT( (ELL-Y)*RS/(A+ELL-Y) ) 

DO 9001=1, N1 
X=(I-1)*DEL 
XPX= XP- X 

IF(XPX.LE.DEL) GO TO 500 
I N ( I , J ) =1 
GO TO 560 

500 IF (XPX .GT.0.0) GO TO 512 

501 IF (XPX. EQ. 0.0) GO TO 515 
IF (J.EQ.l) GO TO 860 

510 GO TO 900 

512 I N ( I ,J)=1 
INC 1 + 1, J ) = 2 
P3 ( J ) =XPX/ DEL 

IF ( J. EQ. 1) GO TO 860 

513 GO TO 560 
515 I N ( I , J ) =2 

GO TO 1000 

560 YPY=ELL-A-X**2/(RS-X**2)-Y 
IF (YPY.GE.DEU GO TO 850 
INC 1+1 » J+l)=2 
P4( I )=YPY/DEL 

850 DIST=XPX*YPY/SQRT(XPX**2+YPY**2) 

852 P4C 1 ) = 1 . 0 

IF (J.EQ.l) GO TO 860 

853 T ( I , J)=TINT*(1.0-EXP(-10.0*DIST*DIST) ) 
GO TO 900 

860 T( I, J )=TINT 
900 CONTINUE 
1000 CONTINUE 
RETURN 
END 



SUBROUTINE TNEW ( N, M, DEL , T I NT , TSURF , T , I N , DELTAT , P4 , P3 , 
* TIM I NT ,F,C ,RHO, ELL,R,S IGMA, FINCR) 

REAL*4 T(20,50) ,P3(50) ,P4(20) , DELTAT (20, 50 ) 

I NTEGE R*2 IN(20,50) 

M1=M+1 

N1=N+1 



SET UP THE TEMPERATURES TO BE USED IN THE MOLECULE 



r ■»'. ^ X 






^ <JU JL» «J 

^ N /j4 ^ 



IF THE TEMPERATURE IS ON THE EDGE CALL TEDGE OTHERWISE 
CALL DERI V 
5 DO 1000 J = 1 » Ml 

10 DO 900 1=1, N1 

IF ( IN ( I , J ) . EQ . 2 ) GO TO 850 

11 IF ( IN ( I , J ) . EQ .0 ) GO TO 860 

12 ZO = T ( I ,J) 

IF ( I . EQ . 1 ) GO TO 15 
Z1=T ( I -1 , J ) 

GO TO 20 
15 Z1=T ( I +1 , J ) 

20 IF (J.EQ.l) GO TO 25 
Z2=T ( I , J — 1 ) 

GO TO 30 
25 Z2 = T I NT 

30 IF ( IN( I ,J + 1).EQ.2) GO TO 50' 

Z4=T( I , J+l) 

GO TO 60 
50 Z4=T ( I + 1 , J+l ) 

60 Z3=T ( I + 1 , J ) 

IF(J.NE.M) GO TO 70 

62 IF(I.EQ.l) GO TO 65 

63 Z4=T ( 1+1, J+l ) 
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GO TO 70 
65 Z4=T ( I , J+l ) 

70 CALL DERI V ( Z0 , Z 1 , Z2 , Z 3 , Z4 , ELL , R ,DE L , TI MI NT , I N , F , 
*C,RH0,P4,P3, DELTAT, I , J , N , A , T , S I GMA ) 

840 GO TO 900 

850 CALL TEDGE(IN,DELTAT,P3»P4,DEL,I,J,T,R,TIMINT,M,F,C, 
*RHC, FI NCR, TINT, SIGMA) 

C KEEP THE BASE AT CONSTANT TEMPERATURE 

GO TO 900 

860 DELTAT ( I , J)=0.0 
900 CONTINUE 
1000 CONTINUE 



• «JL» JLr ^ JU JU JU ou Ju ^ ^ Ur U; JU Uf U/ <JL» U* U/ w U> U> X jUU-UjU/ *jl» gu Ur u> u> u# gu a* u> U# gu gu JU 

^ ■HT 'T' 'T 'T* - T % OT O” Or O^ ‘HP* OP “T* O^ o** O^ rr OT OT OT 'r OT ^r OT ^ -*p> O^ O^ - T" OT 'V* V '1' O^ O^ O* OT OT OT OT O v OT OT OT *T 



INTERPOLATE TO FIND CHANGE AT POINTS WHERE 
BOUNDARY CROSSES Y GRID LINE 

nr nr nr 'r- n" n* nr n" n" nrnrnrn - 'nrnrnrn'' , '|4nrnrnr 44 4r n" n" n- nr n-~ n" 44 nr n" n - * n" 44 n - * •'T' 4r 4 1 ' -v ^ 'l' 'i' t- 4 'i' 'T 4* nr 'T' -r- n~ n" 



DO 2000 J=1,M 
DO 1900 1=2, N1 
IF (I N( I , J) .EQ.0) GO TO 2000 
1700 IF C IN( 1-1, J) .NE.2) GO TO 1900 

1800 DELTAT (I , J ) =0ELT AT ( I , J-l ) +P4 ( 1-1 )* ( DELTAT ( I-1,J)- 
*DELTAT ( I , J-l) ) 

GO TO 2000 
1900 CONTINUE 
2000 CONTINUE 

Jf 0 -» 4 * U« *X> 4 * Uo OL> .t# 4 « <JL> sj- sL U« j, sL> JU J# J# J# U> w *JL» J« J# J# J< J# <JU 4 L» J. -w U# X X X X X X X X X X X X X X X X X X 



CALCULATE THE NEW TEMPERATURE 






jL u* %u U/ u> J# u# u# %l> u> *jju gu gu 
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DO 1010 J=1,M1 
DO 1009 1=1, N1 

1008 T(I,J)=T(I,J)+DELTAT(I,J) 

1009 CONTINUE 

1010 CONTINUE 
1020 CONTINUE 

RETURN 

END 



SUBROUTINE DERIV ( 10 , Z 1, Z2, Z3 , Z4 ,ELL ,R ,DEL , T I MI NT , I N , 
*F,C,RH0,P4,P3, DELTAT, I , J , N , A, T , S IGMA ) 

REAL*4 T( 20,50) ,P3(50) , P4 ( 20 ) , DELT AT ( 20 , 50 ) 

INTEGERS IN ( 20, 50 ) 



X X X X X <J. X X X X X X X X X X X x x x x x X x X X X x X x X X X X x X X X X X X X <JL> X X X X X X X X- x X> x 
X4 v -I- -X nr *t* «•»■* nr nr **» n* •*!' * 1 "* n* *v* *v* nr n'“ n - * 4' 4 s 'r ^ ■■»,-» • ? p 44 -V" n- n~ ■nr n~ *r *r ¥ n* ^ t *r* v *4 t * 1 ' 'r n* v o* t 'f' 'i' 

CALCULATE THE TEMPERATURE DIFFERENCE OVER THE TIME 
INTERVAL 

X X X X X> J< 0- X X X X J - X X X X X X- X X X X X X X X X X X X X X X X X X X X si- X X X X X X sJ- X x x x X X 



IF1J .EQ.l) GO TO 95 
10 N1-=N+1 

Y=( J-l )*DEL 
X= ( I -1 ) *DEL 

I F ( INI 1 + 1, J) .EQ.2) GO TO 40 
35 PP3= 1 . 0 
GO TO 50 
40 PP3 = P3 ( J ) 

50 IF! INI I , J+l) .EQ.2) GO TO 60 
55 PP4 = 1 . 0 
GO TO 70 
60 PP4=P4 ( I ) • 

70 ALPHA=F/ (C*RH0) 
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IF (X.GT.O.O) GO TO 80 
75 X=DE L 

80 Sl= l 2. 0/ ( DELUDED )*(Z3/ ( PP3* ( PP3 +1 . 0 ) J-Z0/PP3+ 

*Zl/l PP 3+1.0 )+Z4/(PP4* (PP4+ 1.0) J-Z0/PP4+Z2/ ( PP4+1.0) ) 
S2=l 1.0/ (X*DEL*(PP 3+1.0) ) ) *( Z3/PP3-PP3*Z 1-Z0/PP3+PP3*Z 

89 ALPHA = F /P.HO 

90 DELTATl I , J ) =ALPHA* ( S1+S2) *TIMI NT 
GO TO 120 

95 DELTAT l I , J)=0.0 
120 RETURN 
END 



SUBROUTINE TE DGE ( I N, DE LT AT , P3 , P4 , DEL , I , J , T , R , T IMI NT , M, 
*F»C»RHO,FINCR, TINT, SIGMA) 

REAL*4 T l 20, 50) ,P3(50) ,P4(20) , DELTAT (20, 50 ) 

INTEGER* 2 IN(20,50) 

M1=M+1 

IF(J.EQ.l) GO TO 900 

C FIND THE X AND Y COORDINATES OF THE POINT 

80 I F ( I.EQ.l) GO TO 550 
40 IF(IN(I-1, J) .NE.2) GO TO 50 
X=( 1-2 )*DEL 
Y=( J-2+P41 1-1) ) *DEL 
45 GO TO 110 

50 IF(J.NE.Ml) GO TO 100 
70 X= ( 1-2 ) *DEL 

Y=( J-2 + P4( 1-1) ) *DE L 
90 GO TO 110 
100 X=( I-2+P3 ( J) )*DEL 
Y=(J-1)*DEL 






CALCULATE THE X AND Y COORDINATES OF THE EDGE POINTS 
ON EITHER SIDE OF THE POINT UNDER CONSIDERATION 






-T- -r- -T" ' 



110 X2=X 

Y3= ( J— 2 ) *DEL 

X3= ( I-2+P3 ( J— 1 ) ) *DEL 

IF(J.NE.Ml) GO TO 106 

Y2= ( J-2+P4 ( 1-1 ) )*DEL 

X3= ( 1-1 ) *DEL 

Yl={ J-2+P3 ( 1-2) ) *DEL 

XI = { 1-3 ) *DEL 

I F ( INC I ,J-1).EQ.2) GO TO 105 

104 Y3=( J-2+P41I ) )*DEL 
GO TO 125 

105 Y3=(J-2)*DEL 

X3= ( I— 2+P3 ( J-l ) ) *DEL 
GO TO 125 

106 IF(IN( I-lt JJ.EQ.2) GO TO 120 
112 IF ( INC I+1,J).EQ.2) GO TO 121 

111 I F ( IN ( I— 1, J+ 1 ) . EQ. 2 ) GO TO 116 

115 Xl=( I-2+P3 ( J+l ) )*DEL 
Yl={ J*DEL) 

Y2=( J-l )*DEL 
GO TO 125 

116 Xl=( 1-2 ) *DEL 
Y1=(J-1+P4( I — 1 ) )*QEL 
Y2 = ( J - 1 ) * D E L 

GO TO 125 

120 XI = { I-3+P31 J))*DEL 
Yl= ( J-l ) *DEL 
Y2=(J— 2+P4( I — 1 ) ) *DEL 
GO TO 125 

121 X3=(I-1)*DEL 

Y3= ( J-2+P4 ( I ) ) *DEL 
Y2 = ( J — 1 ) *D EL 

I F ( IN( I- 1 , J + l) .EQ.2) GO TO 123 
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122 

123 



125 



130 

200 

365 

375 

320 

310 

314 

370 

360 

646 



645 

601 

606 



605 

650 

660 



640 

690 

700 



C 

203 



Xl=( I-2+P3(J+l) ) *DEL 
Y1=J*DEL 
GO TO 125 
Xl=( I-2)*DEL 
Y1=(J— 1+P4( I— 1 ) )*DEL 

# * * Jf if :£ * Jj: •£ ;? :£ * * 'Jf. ^ # # # £ £ ■%. 

CALCULATE THE SLOPE 

sU -Jf %l> U* <X> <JU sU *1# «JL» JU <Xi «l# X X X X X X X X X X X x X X X X X X X X X X X X X X X X X X X X X. X X X X X X X x 



X3X2=X3— X2 
X1X2=X 1-X2 

Vl=Y3-( (X3X2*X3X2) /(X1X2*X1X2) )*Y1 
V2=Y2-( ( X3X2-X3X2) /( XI X2*X1 X2 ) )*Y2 
V3=X3-X2-(X3X2*X3X2)/(X1X2) 

DYDX=( V1-V2) / V3 

X X X X X X X X* X X X X X X X X X X X X X X X X X X X X X X X X X "X X X -X X X -X -X X X X x x X x X X X X X X 
t' ^ ^ ^ ^ ^ ^ ^ ^ <r X ^ ^ ^ ^ ^ n" ^ ^ x ^ X ■¥■ T “T •'r -r 'i' “T *1' 'r -r* *>' 'r -r 'i' 'T' • i r -r -r -r •i' 



FIND THE IMAGE TEMPERATURE 

sU 4f **JU ^ sJU *JL» sJL» si* %L< sL JU %JU «JU %Lr «JU %L* sL sJU JU vL <Ju JU si* -JU <JU sL *1» *JU %L «JU <X <JU «JU \JL» <Ju JL* si* «jU %i> <4# 

-T- -r* -T^ -T^ -T^ -T- -T^ ^ ^ -T^ -T^ >r -T^ -T^ -T^ -T* r^-T^T-n^ -7* v^r 5 ? 



DY=DEL*(SQRT( 1.0/( 1 . 0+DYDX*DYDX ) )) 

DX=DEL* ( SORT (1.0/ ( 1.0+ 1. 0/ (DYDX*DYDX) ) ) ) 

XTERP=( X— DX) — I 1-3) *DEL 

XT=XTERP/DEL 

IF(J.EQ.Ml) GO TO 370 

IF( IN( 1-1, J).EQ.2) GO TO 310 

YT ERP = DEL-DY 

GO TO 360 

IFIDYDX.GT.-2.10) GO TO 370 
YT ERP= P4 < 1-1 )*DEL-DY 
GO TO 360 

YTERP= ( 1 . 0+P4( I— 1 ) ) *DEL-DY 
YT=YTE RP/DEL 

I F { INC I - 1 , J ) .NE. 1) GO TO 645 
GO TO 660 

IN AREAS OF HIGH CURVATURE A DIFFERENT SET OF POINTS 
MUST BE USED 

IF(DYDX.GT .—2.1) GO TO 606 
GO TO 660 
T01=T( 1-2, J— 1 ) 

T1 1=T ( I —1 , J— 1 ) 

T00=T( 1-2, J-2) 

T10=T( 1-1, J— 2) 

TONI =T ( I— 2 , J-3 ) 

IFII.LE.3) GO TO 650 
TN10=T ( 1-3, J-2) 

GO TO 700 
TN10=T( I— 1 , J— 2) 

GO TO 700 
T00=T ( 1-2, J-l) 

T10=T( 1-1, J-l) 

TO 1=T ( 1-2, J) 

TN10=T( 1-3, J-l) 

T11=T( I — 1 , J ) 

IF ( J . EQ.2 ) GO TO 690 
TONI =T ( 1-2, J-2) 

GO TO 700 
T0N1=T INT 

W 1 = T 00* ( 1. 0— XT*XT— YT*YT+XT*YT) 

W2=T 10* ( 0. 5*XT+0.5*XT*XT-XT*YT) 

W3 =T Ml 0* (~0.5*XT+0.5*XT*XT ) 

W4=T01* ( 0. 5*YT+0. 5*YT*YT— XT*YT ) 
W5=T0N1*(-0.5*YT+0.5*YT*YT) 

W6=T 1 1 * ( XT*YT ) 

T I MA GE =1; 1 + W2+ W3+ W4+ W 5+ W6 
SET UP THE EDGE MOLECULES 
Z0=T( I , J) 
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TIMAG2=T IMAGE-2. 0*DEL*SI GMA* (T I MT-T ( I , J) ) 

Z1=TIMAGE 
Z2=T ( I , J-l) 

Z3=T I MAG 2 
204 Z4=T ( I ,J + 1 ) 

PP1=1. 0 
PP3= 1.0 

210 PP2= ( S QRT ( (X2-X3J*(X2— X3) + ( Y2— Y3 J* ( Y 2- Y3 ) ) J/DEL 

211 PP4= ( SORT ( (X1-X2J*(X1-X2)+(Y1-Y2 J*(Y1-Y2)J I /DEL 
IF(J.NE.Ml) GC TO 240 

T I MAG2=T IMAGE-2. 0* DEL* SI GMA*(TINT-T ( I , J) ) 

241 ZO=T ( I ,J) 

Z1=T IMAGE 
• Z3=T I M AG2 

IF( IN( 1-1, Ji.NE.Oi GO TO 349 
345 Z4=T ( I -2 » J ) 

GO TO 350 

349 Z4=T ( I -1 » J ) 

350 IF( IN( I , J-1J .EQ.2J GO TO 207 
Z2=T ( I +1 , J ) 

GO TO 208 

207 Z2=T ( I , J-l ) 

208 PP1=1.0 
PP3= 1. 0 
PP2 = 1 . 0 
PP4=1 • 0 

240 IF( INC 1-1, J) .EQ.2) GO TO 525 
300 IF ( INC 1 + 1, JJ.EQ.2) GO TO 500 
302 GO TO 600 
500 Z2=T(I+1,J) ' 

GO TO 600 
525 Z4=T ( I -1 , J ) 

GO TO 600 
550 ZO=T ( I , J ) 

Z1=T( 1+2, J) 

Z3 = Z1 

Z2=T( I ,J-1 ) 

T I MAG2=T IMAGE-2. 0*DEL*SIGMA* ( TI NT-T ( I , J) ) 

Z4=T IMAG2 
X=l. 0 

590 PP1=1.0 
PP2= PP 1 
PP3 = PP 1 
PP4=PP 1 

600 Sl= ( 2. 0/ ( DELUDED )*( ( 1 . 0/ ( PP3*PP 3+P P3*PP 1) )*( Z3-Z0*( 1. 
*+PP3*Zl/PPl) ) 

S2= ( 2. 0/( DELUDED )*( ( 1.0/( PP4'-PP4+PP4*PP2) )*(Z4-Z0*(1. 
*+PP4*Z2/PP2) ) 

S3=-FI NCR 
ALPHA=F/RHO 

800 DELTAT ( I , J ) =ALPHA* (S 1+S2 ) *T IMI NT+S3 
GO TO 1000 
900 DELTATt I, J)=0.0 
1000 RETURN 
END 
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